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ABSTRACT 
Distributed systems have grown in popularity due to the rapid increase in networking of 
personal computers. A mixture of computers consisting of different architectures can be 
more powerful, reliable, and scalable than a single supercomputer. The problem of 
optimally scheduling jobs on a cluster of heterogeneous machines to minimize the time at 
which the last machine finishes is NP-complete. Nonetheless, the choice of a heuristic 
algorithm greatly affects the speed of solution. This work evaluates a greedy algorithm, 
an A* algorithm, and a simulated annealing algorithm applied to the heterogeneous 
scheduling problem with deadline and dependency constraints. Tradeoffs of speed and 
schedule quality were noted between the algorithms. The greedy algorithm produced 
results quicker than the A* and simulated annealing algorithms, but with a lower 
scheduler quality. Because of these offsetting performance criteria, an analysis was 


conducted to determine which algorithms should be used for which input cases. 
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I. INTRODUCTION 
A distributed system is a collection of heterogeneous computers and processors 
connected via a network [GALLIOO]. This collection works closely together to 


accomplish a common goal. The goal of a distributed operating system is to provide a 





common, consistent global view of the file system, name space, time, security, and access 
to resources. To provide this common view, there are numerous constraints and 
requirements on all participating systems; thus, distributed systems are generally referred 


to as tightly coupled systems. 





Distributed systems have grown in popularity due to the rapid increase in 
personal-computer computing power per dollar. A mixture of computers with different 
architectures can be more powerful, reliable, and scalable than a single big 
supercomputer. A heterogeneous mixture of computers can be more powerful because of 
the sheer quantity of processors in the distributed system; it can be more reliable because 
failure of one computer does not result in failure of the system; and it can be more 


scalable because adding or removing computers from the system is easily accomplished. 





To maximize performance in cluster environments, users can submit programs to 
a resource scheduling and management system that will execute them at the best location, 
retuming the results to the user. Preferably, this management system will execute the 
application in a transparent manner; that is, the user will not know, except perhaps by 
getting better performance, that the application is not executed on the local machine. A 


scheduler for such a heterogeneous environment must be sophisticated. It must 





effectively schedule all of the jobs across the available machines, taking into account 


numerous factors, such as the deadlines of each job, the interdependencies among the 
jobs, and the heterogeneity of both jobs and machines. 

Powerful, reliable, and scalable distributed systems have many applications in our 
current military. Examples of oainent military research in distributed systems are in real- 
time systems like telemedicine, video conferencing, distributed interactive simulations 
such as virtual battlefields, anti-ballistic missile systems, and communication systems. 

The problem of optimally scheduling jobs on heterogeneous machines to 
minimize the time at which the last machine finishes is NP-complete. Nonetheless, the 
choice of a heuristic algorithm greatly affects the speed of solution. This work evaluates 
several approaches for the heterogeneous scheduling problem with deadline and 
dependency constraints. Specifically, three algorithms have been constructed and tested 
for the problem of scheduling jobs on a cluster of heterogeneous machines: a greedy 
algorithm, an A* search algorithm, and a simulated-annealing algorithm. The algorithms 


could be used by the Management System for Heterogeneous Network (MSHN) research 


program currently in progress at the Naval Postgraduate School. 











II. PREVIOUS ATTEMPTS TO SOLVE THIS PROBLEM 


A scheduling problem occurs whenever there is a choice as to the order in which 
tasks can be performed and/or in the assignment of tasks to servers for processing. A 
problem may involve jobs that need to be processed in a manufacturing plant, bank 
customers waiting to be served by tellers, aircraft waiting for landing clearances, or 
program tasks run on a parallel or distributed computer. There is a fundamental similarity 
to scheduling problems regardless of the difference in the nature of the tasks and the 
environment. 

In the era of parallel and distributed computing, scheduling has started to gain the 
attention of many researchers. A program can be viewed as a collection of tasks which 
may run serially or in parallel. The goal of scheduling a set of tasks is to determine an 
assignment of tasks to processing elements and the times at which tasks are executed to 
optimize some performance measures. 

If there are no precedence relations among the tasks forming a program, this 
problem is known as the task allocation problem [ALI94]. Task allocation has been 
studied extensively for the past two decades. The general problem of task scheduling is 
one of the most challenging problems in parallel and distributed computing. It is known 
to be NP-complete even in several restricted cases. Researchers have studied restricted 
forms of the problem by constraining either the task graph or the computer model. These 
cases can be solved efficiently, but are too far from real-world applications. To solve the 


problem in the general case, a number of heuristics have been introduced. These 








heuristics do not guarantee an optimal solution to the problem, but they try to find near- 
optimal solutions most of the time. 

The scheduling investigated in this thesis is static scheduling, the assignment of 
tasks to processors before program execution begins. Static scheduling methods are 
nonpreemptive, meaning that once a job has begun execution, it cannot be interrupted. 
The goal of static scheduling is to minimize the overall execution time of a given set of 


tasks. Work on static scheduling methods by other researchers has attempted to 


[HURSON95}: 
e Predict the program execution behavior at compile time; 


e Group smaller tasks into coarser-grain processes to reduce the communication costs; 
and 


e Allocate processes to processors. 
My work continues with this research, concentrating on the first and third of the bulleted 
items above. 

Static scheduling can be either optimal or suboptimal [HURSON95]. One of the 
most critical shortcomings of static scheduling is that generating optimal schedules has 
been proven, even for it, to be NP-complete except in restricted cases (for example, when 
the execution time of all tasks is the same and only two processors are used). So most of 
the research and development has been focused on suboptimal solutions. These methods 
are classified into exact and heuristic suboptimal approaches. In heuristic methods, the 


solution space is searched in either a depth-first or a breadth-first fashion and the 


algorithm stops when a “good” (or acceptable) solution is reached. 




















The algorithms I present in this thesis illustrate both the exact and heuristic 
approaches. The “greedy” algorithm I implement uses a heuristic approach: A higher 
priority in the scheduling of tasks is given to those that have the smallest estimated 
computation time. The rule of thumb is that by scheduling a task on a machine with a 
smaller than average estimated time, we map the tasks to the machines which execute 
them most efficiently. In contrast, the A* algorithm I implemented uses an exact 
approach: Utilizing both a cost function which measures the actual time expended along 
a path, and an evaluation function which is an estimate of the time from any state to a 
goal, it is possible to find the quickest path to a goal whenever it exists. In order to take 
advantage of both the speed that heuristic approaches provide and the thoroughness that 
exact approaches provide, the simulated annealing algorithm I implemented uses the A* 
algorithm as a base and adds some random heuristics to enhance performance. The 
heuristics will cause the search to explore additional parts of the search graph which may 
result in finding a goal to the path quicker than the A* algorithm did. 

Static scheduling suffers from a wide number of additional problems, most 
notably [HURSON95]: 

e The lack of efficient and accurate methods for estimating task execution times and 
communication delays can cause unpredictable performance degradations. Compile- 
time estimation of the execution time of a program’s tasks is often difficult due to 
conditions or iteration counts that are unknown before execution. 

e Existing task scheduling methods often ignore that data can be geographically 
distnbuted. This omission causes performance degradations due to run-time 


communication delays for accessing data at remote sites. 


e Static scheduling schemes need a tool to provide a performance profile of the 
predicted execution of the scheduled program on a given architecture. The user can 


then utilize this tool to improve the schedule or to experiment with a different 
architecture. 


Although my thesis does not address these problems, they are areas researchers are 
currently studying. Several students and faculty performing research on the iceaaea 
System for Heterogeneous Network (MSHN) have done research on these topics; they 
are critical to the success of the scheduling component of MSHN. Armstrong, Hensgen, 
and Kidd investigated estimation of task execution times, and discussed the relative 
performance of various mapping algorithms [ARMSTRONG98]. Bhat, Raghavendra, 
and Prasanna investigated various aspects of communication delays and their effect on 
performance in a distributed system [BHAT99]. Alhusani, Prasanna, and Raghavendra 
investigated a unified resource scheduling framework for heterogeneous computing 


environments [ALHUSAINI99]. Similar work was done with MSHN 


[FREUND98],[KIDD96]. 








il. DESCRIPTION OF THE APPLICATION 


A. SMARTNET 

The research which led to the current MSHN project originated with the SmartNet 
project. SmartNet was a scheduling framework for heterogeneous computing. Several 
research groups contributed to its development: NCCOSC RDTE Division in San Diego, 
SAIC, the University of Cincinnati, and the Naval Postgraduate School. SmartNet was 
different than other scheduling frameworks. Most frameworks use policies that employ 
Opportunistic Load Balancing (OLB). SmartNet, however, took into account not only the 
loads on various machines, but also, the performance of the application on various 
platforms (heterogeneity) [MSHN0O0]. 

The scheduling engine of SmartNet employed several scheduling algorithms. 
Two algorithms, “gensked” and “metasked”, took into account scheduling of jobs with 
precedence constraints. Gensked and metasked belong to a class of “generational” 
scheduling algorithms that schedule the jobs starting at the root of the directed acyclic 
graph of task precedences and schedule a task after all its predecessor tasks are scheduled. 
They both use static algorithms. 
B. MSHN 

The MSHN project is an extension of the SmartNet research. The goal of MSHN 
is to develop a distributed system which will have naval shipboard applications. It is a 
part of the DARPA/ITO QUORUM program which is developing iechinolonies that allow 


end users to achieve predictable and controllable end-to-end quality of service (QoS) for 





critical defense computing needs in a global heterogeneous distributed computing 
environment. The MSHN team is a collaboration between DoD (Naval Postgraduate 
School), academia (Naval Postgraduate School, University of Southern California, 
Purdue University), and industry (NOEMIX). The MSHN architecture will permit 
execution of multiple, disparate tasks that use a set of shared, heterogenous resources ina 
way that maximizes a collection of application-specific quality of service (QoS) 
measures. Factors influencing QoS for MSHN include security, deadlines, priority, 
adaptability (preferences for different versions), resource availability, and external users. 

MSHN’s components will perform several key functions. These include (1) 
monitoring general resource availability, (2) transparent sensing of resource requirements 
of applications, (3) on-line measurement of resource and system state, (4) mapping tasks 
and subtasks onto a heterogeneous suite of machines in a way that exploits heterogeneity, 
(5) adaptation of jobs to variations in the availability of resources, and (6) meeting QoS 
requirements including real-time deadlines, fault tolerance, security, and priorities. 
Because MSHN will be operating in an environment where usages and loads are not 
predictable, uncertainty will be quantified and accounted for in all of MSHN’s 
components. 

MSHN’s architecture, shown in Figure 1, consists of (1) a client library that wraps 
system calls (including process-start requests), (2) a daemon (mshnd), (3) the Scheduling 
Advisor, (4) the Resource Requirements Database (RRD), and (5) the Resource Status 


Server (RSS). Users may submit jobs at any of the client sites. Process-start requests are 


sent to the Scheduling Advisor. The Scheduling Advisor queries the RRD for historical 











ms ) 
Advisor 
ystems 


apl. |} 
Scheduling 
Operating Sy 
Figure 1. MSHN Architecture. 
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and other information on the job, and also the RSS for information on status of 
resources. The Scheduling Advisor then determines the resources that are to be allocated 
for the job. The job is executed on these resources. Finally, the RRD and the RSS are 
updated. MSHN may use CORBA in its architecture, as well as integrating the above 
scheduling functions into an ORB (Object Request Broker) [MSHNOQ]. 

Scheduling is important to the MSHN project since it particularly addressed 
military applications and these generate scheduling problems. The detection and tracking 
of aircraft is a problem which needs efficient scheduling to ensure reliable performance of 
the system. Resources can’t be solely dedicated to detection, or solely dedicated to 
tracking one platform; scheduling must allow adequate allocation of resources such that 
each part of the air traffic control problem is satisfied. Or consider a naval ship which is 
required to detect, track, and target incoming missiles. Efficient scheduling is necessary 


to adequately allocate resources such that each part of the naval ship’s air defense mission 


is satisfied. 
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IV. DESCRIPTION OF THE PROGRAM 


A. NOTATION 

A task is a program or a set of programs to be executed upon a single request from 
auser. A task can consist of multiple jobs that can be distributed to multiple machines. 
Precedence constraints among jobs can be represented by ordered pairs of the form (Ji > 
Jj), which means that job Jj must not start executing before job Ji completes. Therefore, 
the precedence constraints define a partial order on the set of jobs that constitute a task. 
Hence, we can represent the tasks and precedence constraints by a Directed Acyclic 
Graph (DAG), which need not be a tree. We define the set of jobs that must complete 
before a job j can begin to be the Pred(j). Similarly, we define j to be a ready job if all of 
the jobs in Pred(j) have completed. 

Figure 2 illustrates the above terms. A DAG representing a task consisting of 15 
jobs is shown. The precedence constraints are shown by directed edges. Initially, Jobs 
JO, J1, and J2 are ready jobs. Jobs J3, J4, and J5 cannot be started until JO completes. 
Similarly, Job J9 has to wait for both J3 and J4 to complete, before it becomes ready. 

We now define an ETC matrix. ETC(j,m) is the expected time to complete job j 
on machine m, in a sterile environment. A sterile environment is one where no other jobs 
for the task use any other machines or the network connecting them. 

A schedule is an assignment of the jobs to the available machines such that 
(1) the precedence constraints are not violated; 
(2) ajob runs on exactly one machine; 


(3) a machine runs no more than one job at the same time; and 
(4) deadlines are met. 


I] 


eady Jobs 


Precedence constraint 


Pe oe 


Figure 2. A Directed Acyclic Graph (DAG) Representing a Task. 














The schedule will specify, for all jobs: 


(1) The machine on which the job should be executed; 
(2) The expected start time of the job on the selected machine; and 
(3) The expected completion time of the job on the selected machine. 


B. THE SCHEDULING PROBLEM 

The classic job scheduling problem is to determine an optimal schedule for a set of jobs 
and a set of machines, given precedence and deadline constraints, if any. For example, 
consider the DAG and the corresponding ETC matrix shown in Figure 3. Two different 
schedules are also shown in Figure 3 which satisfy this DAG. The DAG shows that job 2 
and job 3 can start only after job 1 finishes. Similarly, job 4 can start only after both job 2 
and job 3 finish, job 5 can start only after job 3 finishes, and job 6 can start only after 
both job 4 and job 5 finish. Its corresponding ETC matrix lists the ETC values of the jobs 
1 through 6, for four different machines. For example, job 1 is expected to complete after 
36 time units if it is executed on machine C, and job 3 after 41 time units if it is executed 
on machine A. 

Figure 4 shows two example schedules for the DAG of Figure 3. Each row of the 
schedule depicts the execution of the machine corresponding to its label. Shaded boxes 
represent jobs. The job numbers are labeled. Empty boxes represent machine idle time. 
Schedule A has a length of 81 time units, and Schedule B, 90 time units. Schedule A is 
an optimal schedule. The best possible completion time will be bounded below by the 
time of the longest subsequence that starts at the root node of the DAG and ends at a leaf 


node. There are three sequences 1-2-4-6, 1-3-4-6, and 1-3-5-6 in this example DAG, | 


13 


having lengths 65, 73, and 81, respectively. The length of the longest path is 81. Hence, 

schedule A is an optimal schedule. 

We note that machines can be idle in an optimal schedule, since there are 
precedence constraints among the jobs. We also note that if multiple tasks are submitted 
to the system, there will be multiple DAGs. However, throughout this thesis, we assume 
without loss of generality that there is only one DAG representing all the tasks since we 
can construct a master DAG which consists of a dummy root node and has the top level 
nodes of the DAGs representing the tasks as its children. See Figure 5. 

Each job is assumed to have a deadline, a time by which the job should be 
completed. The algorithms implemented use a “soft” deadline: If a job does not 
complete by its deadline, the system does not fail. However, each algorithm uses a cost 
function which penalizes missed deadlines. The more time from the deadline of the job 
to the start of the job, the higher the cost of the job. This relationship is linear. 

It is useful to identify some additional assumptions made here. 

e These algorithms were tested in a linear eterosencous environment where the ratio 
of the processing powers of any two machines is a constant. If two machines M1 and 
M2 are linearly heterogeneous, that means that one job running p times faster in 
machine M1 than in machine M2 implies that all jobs run p times faster in machine 


M1 than in machine M2. 


e These algorithms are for deterministic scheduling where it is assumed that complete 
information about the job to be scheduled is available prior to execution. 
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Figure 3. Directed Acyclic Graph (DAG) in a Heterogeneous Environment. 
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Schedule A completed at time 81 





Figure 4. Multiple Schedules Satisfying the Same Directed Acyclic Graph (DAG). 
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Figure 5. A Master Directed Acyclic Graph (DAG) representing Multiple Tasks. 
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These algorithms are static scheduling algorithms. Scheduling algorithms can be 
classified based on the time at which scheduling decisions are made. In static 
scheduling algorithms, information regarding the entire task graph must be estimated 
prior to execution and all jobs are assigned to machines before any job starts 
execution. 


These algorithms use non-preemptive scheduling, meaning that once a job has begun 
execution, it cannot be interrupted. 


GREEDY ALGORITHM 


A “greedy” algorithm was implemented to produce a schedule by a limited number of 


- uncomplicated heuristics. Data structures used in the implementation of the algorithm 


were. 


ETC[][]: A two-dimensional integer array is used to store the ETC values. As 
discussed, ETC(j,m) is the expected time to complete job j on machine m, assuming a 
sterile environment. 


DEADLINES[]: A one-dimensional integer array is used to store the deadlines for 
each job. 


DEPENDENCIES [][{]: A two-dimensional integer array is used to represent 
precedence constraints among jobs. A “1” entered in the array represents precedence 
constraints between the i job and the j job — specifically, that job Jj must not start 
executing before job Ji completes. We say then that Ji is a parent of Jj, and Jj is a 
child of Ji. 


CHECKSUM[]: A one-dimensional integer array is used to represent the current 
ability for a job to be scheduled without violating precedence constraints. Each index 
in the array corresponds to a job. The value stored represents the number of jobs 
currently required to be completed prior to the execution of that job. If the value is | 


zero, then the job can be scheduled without violating any precedence constraints. 


Initial values are obtained by counting the nonzero values from the corresponding 
column of the dependency array. As jobs are scheduled, the array is updated to reflect 
current precedence constraints. 


MACHINE_ETC[]: A one-dimensional integer array is used to represent the stop time 
of each machine, the time at which the machine will complete the job currently 
scheduled on it. 
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TASK_STATUS[]: A one-dimensional boolean array is used to represent the status of 
jobs scheduled, whether scheduled or not. 


SCHEDULE_SUMMARY{][{]: A two-dimensional array is used to store the results 
of the scheduling. The rows represent the jobs scheduled, and four columns give the 
job number, the machine that job was scheduled on, the start time of the job, and the 
completion time of the job. 


The following steps are followed by the greedy algorithm. 


Have all tasks been scheduled? If NO, continue with scheduling. 


Look in the CHECKSUM array to determine all ready jobs, jobs that have no 
precedence constraints and can be scheduled. 


Look in the TASK_STATUS array to determine which of these jobs have not been 
scheduled. 


Find the minimum ETC value for each job which does not violate the corresponding 
job deadline in the DEADLINE array. 


Select the minimum of these values over the jobs and identify the corresponding job 
and machine. Schedule this job on this machine. 


Update the TASK_STATUS array to indicate the selected job has been scheduled. 


Update the CHECKSUM array if the job scheduled was a parent of any other job 
remaining to be scheduled. 


Repeat the above choices until all jobs have been scheduled. 

The SCHEDULING_SUMMARY array will hold the final schedule. 
A* ALGORITHM 
1. General Description 


The greedy algorithm discussed above is not guaranteed or designed to find the 


shortest schedule. This is desirable for our scheduling problem, as the main justification 


for distributed processing is speed. A heuristic that finds the minimum-cost path to a 
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goal whenever it exists is said to be admissible [LUGER89]. One example of an 
admissible heuristic is that used by a restricted kind of breadth-first search, that explores 
in a level-by-level fashion. Since it looks at every state at level n of the graph before 
considering any state at the level n+ 1, any goal nodes are found along the minimum-cost 
path provided that costs do not vary over branches. 

2. Cost Functions 

A cost function measures the work expended along a path between states. In our 
scheduling problem, the cost from the starting state (no jobs scheduled) to its descendants 
can be measured as total computing time. The cost of scheduling a specific job on a 
specific machine is simply the time it takes that job to execute on that machine. Then the 
cost of a state where multiple jobs are scheduled on multiple machines is the maximum of 
the individual machine completion times, or the time the last job completes. 

3. Evaluation Functions 

An evaluation function is an estimate of the cost from any state to a goal. In our 
scheduling problem, the cost from any state (where some jobs are scheduled) to the goal 
(where all jobs are scheduled) needs to incorporate the ETCs of the unscheduled jobs and 
the processing capabilities of the machines in the network. The best that could be 
accomplished at any given state would be to run all E the unscheduled jobs on the 
available machines with no idle time on any machine. To perform this calculation, we 
need the total estimated run time of all unscheduled jobs, and the total number of 
processors in the network. The estimated run time of an arbitrary unscheduled job is the 


mean of the estimated time to compute values over all machines for that job. The total 
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estimated run time of all unscheduled jobs is the sum of the individual jobs estimated run 
times. The total number of processors in the network is known prior to scheduling. 
Therefore, to model all of the unscheduled jobs running on available machines with no 
idle time on any machine, we divide the total estimated run time of all unscheduled jobs 
by the total number of processors in the network. That total time could be used as a lower 
bound on the cost from any state to the goal. A* algorithms with a lower-bound 
evaluation function are admissible because they always find the optimal solution path as 
the first path to a goal. This is an important attribute for our scheduling problem, as we 
can use this as we compare results from A* to results of other algorithms. 
4, Data Structures 
Data structures used in the implementation of the algorithm are: 


e ETC{[][]: A two-dimensional integer array is used to store the ETC values as with the 
greedy algorithm. 


e DEADLINES[]: A one-dimensional integer array is used to store the deadlines for 
each job as with the greedy algorithm. 


e DEPENDENCIES[][J: A two-dimensional integer array is used to represent 
precedence constraints among jobs as with the greedy algorithm. 


e agenda: A list representing all states visited that have not been evaluated for successor 
states or evaluated to be the goal state. 


e usedagenda: A list representing all states visited that have been evaluated for 
successor states and have been determined not to be the goal state. 


E. SIMULATED ANNEALING ALGORITHM 
The simulated annealing algorithm modifies the A* algorithm. In an A* search, 


the state chosen to be visited is the state with the lowest total cost plus evaluation 


Z1 


function. In simulated annealing, that state has the highest probability of being chosen, 
but the other states are also possible with a probability. For example, if the first state in 
the sorted agenda had a total cost plus evaluation function of 20, and the second state had 
21, the desirability of either state is virtually the same. However, if the first was 20, and 
the second was 50, the first should be definitely more desirable. So the probability of 
choosing an item should depend on its total cost plus evaluation function relative to the 
first item in the agenda. The reason simulated annealing might perform better 1s that an 
A* search might start on a part of the search graph that looks promising, but is not a good 
path to the goal at all, and we don’t know that until we wasted valuable computation 
time. With simulated annealing, occasionally we explore other parts of the search graph. 

As the search continues deeper into the search graph, we are more likely to have 
found the start of the best salle and need to decrease the probability of choosing a state 
other than the first state. Thus, as time progresses, simulated annealing becomes ini 
like A* search. To implement simulated annealing, the likelihood of choosing a state in 
the sorted agenda can be proportional to the ratio of the total cost plus evaluation function 
of the best state on the agenda to the total cost plus evaluation function of the given state. 
For example, if the first state had a total cost plus evaluation function of 20, and the 
second state had a cost plus evaluation function of 40, we would be half as likely to select 
the second state than the first state. Likelihoods can be converted to probabilities by 
dividing by their total. 

As we progress in our search, we want the likelihood of choosing states other then 


the first state in the sorted agenda to decrease. We do this by dividing the likelihood for 
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each item except the first item by the time elapsed in the search. The method used to 
measure the time elapsed in the search is to count the number of steps we have progressed 
in the search. As the search progresses, this will result in a lower and lower likelihood of 
choosing a state other than the first state in the sorted agenda. 
F. TEST SETUP 
Because the algorithms implemented in this thesis were written in different languages, a 
program for generating test cases for input was necessary. The program was written in 
Java and outputs the following information to a text file which can be read by the 
algorithms: 
e Number of jobs to be scheduled 
e Number of machines in the network 
e ETC matrix 
e Deadlines for each job 
e Precedence constraints among jobs 

Three levels of precedence constraints were modeled: low level, medium level, and 
high level. A low level of precedence constraints is achieved by filling the entire 
dependency array with zeros. A high level of precedence constraints is achieved by 
filling the portion of the dependency array below and to the left of the arrays diagonal 
with ones, and filling the rest of the array with zeros. A medium level of precedence 
constraints is achieved by filling the portion of the dependency array below and to the left 


of the arrays diagonal with some ones and some zeros. 
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Professor Neil Rowe provided the framework of the A* algorithm written in LISP. 


Utility functions were written by myself to interact with it. 
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Vv. DISCUSSION OF RESULTS 


The greedy algorithm, A* search algorithm, and the simulated annealing 
algorithm were run to produce schedules. The greedy algorithm was written in Java, and 
the A* and simulated annealing algorithms were written in LISP. We varied the 


parameters listed for the test cases in the last section. The results of the scheduling give: 


e The machine on which the job should be executed 

e The expected start time of the job on the selected machine 

e The expected completion time of the job on the selected machine 
e The time the last job completes over all machines 


e The overhead of the program (how long it took the program to compute the schedule) 


Table 1 shows a summary of these test runs. Two criteria can compare algorithms: the 
time the last job completes, and the program overhead. With relatively simple scheduling 
problems, we see that the schedules produced by the three algorithms provide similarly 
good schedules. As the difficulty of the scheduling problem increases (the number of 


jobs needing to be scheduled increases), the A* search algorithm and the simulated 


annealing algorithm produce better schedules than the greedy algorithm. This was 


expected, as the greedy algorithm is not making necessarily wise choices. But as the 
difficulty of scheduling increases, the computation time increases faster for A* and 
simulated annealing than for the greedy algorithm. Although Table 1 shows no increase 


in algorithm overhead for the greedy algorithm as the difficulty of the scheduling problem 
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increases from 10 jobs to 50 jobs to 100 jobs, an increase was seen with much more 
difficult scheduling problems. 

Another point of interest is comparing A* and simulated annealing. Note that for 
the scheduling problems with high and medium-level precedence constraints, the 
simulated annealing algorithm did not find a solution faster than the A* algorithm. 
However, with no precedence constraints the simulated annealing algorithm did find a 
solution faster than the A* algorithm. The original intent of modifying the A* algorithm 
to a simulated annealing algorithm was to produce an algorithm with behaviors similar to 
an A* search, but which would find solutions in complex search problems faster than A*. 
Whether this has been accomplished is inconclusive at this point. A wider variety of test 


runs or different input data may be required to resolve this question. 
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Table 1. Results of scheduling using a greedy algorithm, A* search algorithm, and 
simulated annealing algorithm with various level of precedence constraints among jobs. 
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VI. CONCLUSION 


The implementation of a variety of algorithms with different heuristics and 
different properties is useful in scheduling. Simple algorithms produce results quicker 
than more complicated algorithms, but the schedules they produce are usually not as good 
as the schedules produced by the more complicated algorithms. Because of these 
offsetting performance criteria, a detailed analysis using a variety of input cases must be 
performed to determine which algorithms should be used for which input cases. The | 
information gained by such an analysis would be useful in the implementation of a 
scheduler. The scheduler could analyze the makeup of a scheduling problem, compare 
that to various test cases, and choose the algorithm to do the scheduling that will perform 


best under initial conditions of the scheduling problem. 
Future work which could extend this thesis could include: 


e Implementing other scheduling algorithms and comparing performance 
e Modeling the ETC values to known or expected inputs of a specific problem 
e Modeling the deadline values to known or expected deadlines of a specific problem 


e Modeling the precedence constraints to known or expected precedence constraints of 
a specific problem 


e Implementation of a GUI-based simulation tool giving the tester a user-friendly 
graphical way to alter input conditions and view scheduling results 
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APPENDIX As 


public class scheduler implements setupConstants { 


static 
static 
static 
static 
static 


static 
static 
static 
static 
static 
static 
static 
static 





TEXT OF THE PROGRAM 


int NUM_TASKS 

int NUM_MACHINES 
long START_TIME 
long FINISH_TIME 
long COMPUTE_TIME 


| | 
“eo 


4 ~e 


“=e 


ane eo > ee ee) 


“es 


Bred tl? 

DEADLINE []; 
DEPENDENCIES [][{]; 

int CHECKSUM[]; . 
int SCHEDULE_SUMMARY([][]; 
String TASK_ORDER[] []; 
boolean TASK_STATUS[]; 
boolean MACHINE STATUS[]; 


a Ie lh es 
Int 
ant 


generateMatrices generate; 
generateTestFiles test; 
generateResultFiles result; 
displayMatrices display; 
search_Greedy greedy; 


[ [BRR ERRERREKRKEREREKEKERKREKEEKKEEKERES 


//ARGUMENTED CONSTRUCTOR 


L [BRR RERKRKKKRKKEKKEKEKREKREKKEKKEKKK KKK 


public scheduler({ int numtasks, 


{ 


SCHEDULE_SUMMARY 


NUM_TASKS 
NUM_MACHINES 


int nummachines, String filename } 


numtasks; 
nummachines; 


new int [NUM_TASKS] [4]; 


generate = new generateMatrices (NUM_TASKS, NUM_MACHINES) ; 

display = new displayMatrices (NUM_TASKS,NUM_MACHINES) ; 

greedy = new search Greedy (NUM_TASKS,NUM_MACHINES) ; 

test = new generateTestFiles (NUM_TASKS,NUM_MACHINES, filename) ; 
result = new generateResultFiles (NUM_TASKS,NUM_MACHINES, filename) ; 
ETC = generate.generateETC (); 

DEADLINE = generate.generateDeadlines(); 

DEPENDENCIES = generate.generateDependencies (dependency_type) ; 
CHECKSUM = generate.generateCheckSum(DEPENDENCIES) ; 

TASK ORDER = generate.initializeTaskOrder (); 

TASK STATUS = generate.initializeTaskStatus(); 

MACHINE_STATUS = generate.initializeMachineStatus(); 


display.displayInitialization (ETC, DEADLINE, DEPENDENCIES, CHECKSUM) ; 
test.createTestFile (ETC, DEPENDENCIES, DEADLINE) ; 


}//END CONSTRUCTOR 


[ [RRR KRKKEEKRKEEKKKKKKKEKKKKKKK KK 


//SCHEDULE THE META-TASK 


[LL RRR RRR RKKEKKEKKKEKEKKKEKKEKKKKKKKK 


public void scheduleMetaTask() 


{ 


3] 


START_TIME = System.currentTimeMillis(); 

SCHEDULE_SUMMARY = 
greedy .search (ETC, TASK_STATUS , MACHINE_STATUS , TASK_ORDER, CHECKSUM, DEPENDE 
NCIES, SCHEDULE_SUMMARY) ; 

FINISH TIME = System.currentTimeMillis(); 

result.createSummaryFile (SCHEDULE_SUMMARY, START_TIME, FINISH_TIME) ; 

}//END METHOD "schedule" 


LL RRRKRERRRKKRERKEKEKREREREKEERERERERER 


//DRIVER 


[ [RRR RK KK KKH KE KEKE KKKEKREKKKEKKE 


public static void main(String args[]) 


{ 
if (args.length < 3 || args.length > 3){ 
System.err.println("USAGE: java TwoDArray <l>#tasks <2>#machines 
<3>filename") ; 


System.exit(1); 
} 


scheduler ETC = new scheduler(Integer.parseiInt(args[0]), 
Integer.parseInt(args[1]), 
args[2]); 
ETC .scheduleMetaTask(); 
}// END METHOD "main" 


}// END CLASS "scheduler" 
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import java.util.Random; 
public class generateMatrices implements setupConstants{ 


public int NUM_TASKS; 
public int NUM_MACHINES; 


//use seed value to allow repeatable experiments 
long seedValue = 0; 
Random r = new Random( seedValue ); 


[LL RRR RRKRKRKEKKKKEKKERKKKRKEKKKKEKE RK KKK 


//ARGUMENTED CONSTRUCTOR 


[ [RRR KKK KKK KKK KEKE KEK ERK KKK 


public generateMatrices(int numtasks, int nummachines) 
{ 

NUM_TASKS = numtasks; 

NUM_MACHINES = nummachines; 
}//END CONSTRUCTOR 


[LL RER RRR RR KRRKKRKEKKEKRKEKKEKEKKKRK KK KK KKK 


//GENERATE THE ETC MATRIX 


[LER RRKRREKRKR KEKE KEK KKK KK KKK KKK KEKE 


public int[][] generateETC() 
{ 


int ETC[][] = new int [NUM_TASKS] [NUM_MACHINES] ; 

for(int task = 0; task < NUM_TASKS; task++) 

for(int machine = 0; machine < NUM_MACHINES; machine++) 
ETC[task] [machine] = Math.abs(r.nextInt()) % ETC_RANGE + 1; 

} 


return: ETC; 
}//END METHOD "generateETC" 


[LL RRRRRKRKRRKEKKEKKEKKEKKKKKKKKEKKEKKKKEE 


//GENERATE DEADLINES 


LL RRR RRR KEKE KEKEKEKEKKEKKKKKKKKEK 


public int{] generateDeadlines() 


{ 
int DEADLINE[] = new int [NUM_TASKS]; 


for(int task = 0; task < NUM_TASKS; task++) 
{ 


} 
return DEADLINE; 
}//END METHOD "generateDeadlines" 


DEADLINE[task] = Math.abs(r.nextInt()) % ETC_RANGE + 1; 


JL RRRRKRKRREKKKRREKEKERKEKKRKREKEKKEKKEKEEE 


//GENERATE DEPENDENCIES 


[LL RR RRR RKRKRKKKKKKKKKKK KKK KKK KKK KEK 


public int[{][] generateDependencies( int dependency_type ) 


{ 
int DEPENDENCIES[][] = new int [NUM_TASKS] [NUM_TASKS]; 
double random_dependency; 


for(int first_task = 0; first_task < NUM_TASKS; first_task++) 
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{ 
for(int second_task = 0; second_task < NUM_TASKS; second_task++) 


{ 
//set diagonal of dependency matrix to zeros 
//--> a task can’t be dependent on itself 
//also set lower half of dependency matrix to zeros 
//--> eliminate the possibility for cyclical dependencies 
if( first_task <= second_task ) 
{ 
DEPENDENCIES [first_task] [second_task] = 0; 


} 


//full dependencies 
else if(dependency_type == FULL_DEPENDENCY) 
{ : 

DEPENDENCIES [first_task] [second_task] = 1; 


} 


//no dependencies 
else if (dependency_type == NO_DEPENDENCY) 
{ 
DEPENDENCIES [first_task] [second_task] = 0; 
} 


//randomly generated dependencies 
else if(dependency_type == RANDOM_DEPENDENCY) 
{ 

random_dependency = Math.random(); 

if (random_dependency < 0.5) 


DEPENDENCIES [first_task] [second_task] = 0; 
else 
DEPENDENCIES [first_task] [second_task] = 1; 


} 
} 


} 
return DEPENDENCIES; 


}//END METHOD "generateDependencies" 


/* 


[ [RRR ERE RRREKREEERERE REESE RERRRER EEE 


//TRANSITIVE CLOSURE OF DEPENDENCIES 


[ [RRR EEERERERERERERERREE ER ERRERE REE E 
public int[][] transitiveClosureOfDependencies( int DEPENDENCIES[][] ) 
{ 
// USE TRANSITIVE CLOSURE ALGORITHM HERE 
}//END METHOD "transitiveClosureOfDependencies" 


oy 
LL ERERERESERRERRERER EERE REE EE REAR 


//CHECK-SUM OF DEPENDENCY MATRIX 


ff FERRER ERRREREREREA EERE AEE BARR 


public int{] generateCheckSum( int DEPENDENCIES[][] ) 


{ 
int CHECKSUM[] =’ new int [NUM_TASKS] ; 


for(int first_task = 0; first_task < NUM_TASKS; first_task++) 
{ 


for(int second_task = 0; second_task < NUM_TASKS; second_task++) 
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{ 
i1£({ DEPENDENCIES [first_task] [second_task] 


CHECKSUM [second_task] ++; 


I 
Hl 
= 


} 
} 
return CHECKSUM; 
}//END METHOD "generateCheckSum" 


[LER RRR KKK KREKRKEKKEKKRKEEKKEKRKEKEE 


//INITIALIZE TASK ORDER 
[BRI RI KR FR IO IR II RI kk 


public String[{][] initializeTaskOrder () 


{ 
String TASK_ORDER[][] = new String[{NUM_TASKS] [NUM_MACHINES] ; 


for(int task = 0; task < NUM_TASKS; task++) 
{ 
for(int machine = 0; machine < NUM_MACHINES; machine++) 
{ 
TASK ORDER[task] [machine] = "NO"; 
} 
j 
return TASK ORDER; 
}//END METHOD "initializeTaskOrder" 


[ [RRR RRKREKEKRKEKKKEKRKKKKKKKKEKKKKEKEKE 


//INITIALIZE TASK STATUS 


[ [RRR RRR RRR KEK KEEKREKKEKEKKEKEKRKEKKKRKES 
public boolean[] initializeTaskStatus () 


{ 
boolean TASK_STATUS[] = new boolean[NUM_TASKS]; 


for(int task = 0; task < NUM_TASKS; task++) 
{ 


} 
return TASK_STATUS; 
}//END METHOD "initializeTaskStatus" 


TASK _ STATUS[task] = false; 


[ [RRR RRR KRRKKE KEKE KRREEKKKEEERKEKRE KR 


//INITIALIZE MACHINE STATUS 


[ [RRR RRRKRRKEKKRKRKREREKKEKEKEKKKEKKKKE 
public boolean[] initializeMachineStatus() 


{ 
boolean MACHINE_STATUS[] = new boolean [NUM_MACHINES]; 


for(int machine = 0; machine < NUM_MACHINES; machine++) 
{ 
MACHINE STATUS [machine] = false; 
return MACHINE_STATUS; 
}//END METHOD "initializeMachineStatus" 


[LL RRR RRR RKRKEKRKKKKEKKKKRKREKKEKKKEKRKE 


//INITIALIZE MACHINE ETC 


[ [RRR RRKEKKKRKEKK KR KK KKKKKEKKEKEKKEEE 
public int{] initializeMachineETC () 


{ 
int MACHINE_ETC[] = new int [NUM_MACHINES]; 
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for(int machine = 0; machine < NUM_MACHINES; machine++) 
{ 
MACHINE ETC[machine] = 0; 
} 
return MACHINE_ETC; 
}//END METHOD “initializeMachineETC" 


}//END CLASS "generateMatrices" | : 
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//creates and writes to a file using printstream object 


import java.io.*; 
import java.awt.*; 


public class generateTestFiles 
{ 

public PrintStream ps; 

IOS os = new I0S(); 

public int NUM_TASKS; 

public int NUM_MACHINES; 


public generateTestFiles(int numtasks, int 
filename) 
{ 
NUM_TASKS numtasks; 


NUM_MACHINES nummachines; 


String file = new String(filename) ; 


bry 4 
ps = new PrintStream ( 
new FileOutputStream(filename) ); 
} 
catch( IOException e ) { 


nummachines, String 


System.err.println("File not opened properly\n" + 


e.toString()); 
System.exit(1); 
} 
} 


public void addTestSetup() 
{ 
ps.printin(NUM_TASKS) ; 
ps.printin(NUM_MACHINES) ; 
pS<printin() ; 


} 
public void addETC( int ETC[][] ) 
for(int i = 0; i < NUM_TASKS; i++) 
for(int j = 0; 3 < NUM_MACHINES; j++) 
os.setJustify(IOS.LEFT) ; 
os.setwidth(5); 
os.print(ps,ETC[i] [j]); 
ere 
ne 


public void addDependencies( int DEPENDENCIES [] [] 
for(int i = 0; 1 < NUM_TASKS; i++) 
for(int j = 0; j < NUM_TASKS; j++) 
os.setJustify(IOS.LEFT) ; 
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os.setWidth(5) ; 
os.print(ps,DEPENDENCIES[i][j]); 
} 
ps. printin(); 
} 
ps.println(); 
} | 


public void addDeadline( int DEADLINE[] ) 
{ 
for(int i = 0; i < NUM_TASKS; i++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth(5) ; 
os.print(ps,DEADLINE[i]); 
DS.peinelal) 
} 
} 


public void createTestFile( int ETC[][], 
int DEPENDENCIES[][], 


int DEADLINE[] ) 


System.out.printin("Generating Test File...."); 
System.out.printin(); 

addTestSetup(); 

addETC (ETC) ; 

addDependencies (DEPENDENCIES) ; 

addDeadline (DEADLINE) ; 
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public class search_Greedy implements setupConstants { 


// *** Performance Metric *** 

// Simply review the estimated time to complete (ETC) matrix 
// and choose the task on the machine which will run fastest 
// *** BRilgorithm Description *** 

// 1. Find the minimum ETC for each task -- 


fy take into consideration also time which each machine is available 
// 2. Choose the minimum of the minimums (i.e., the min-min 

Ti algorithm), and schedule that task 

// 3. Update the task status and machine status to indicate 

a which tasks are scheduled, and which machines are running 


// 4. Do steps 1 thru 3 until all tasks are scheduled 


public int NUM_TASKS; 

public int NUM_MACHINES; 
int[] MACHINE STOP _TIME; 
displayMatrices display; 


[ [RRR RRR REKRKEKEEKKEKRKEKEKR KKK KKK KKK KEKE 


//ARGUMENTED CONSTRUCTOR 

[LL RRR RRR RRR KKKKKR KKK KKK KKK KKK EK 

public search_Greedy(int numtasks, int nummachines) 

{ 
NUM_TASKS 
NUM_MACHINES 
MACHINE_STOP_TIME 


numtasks; 
nummachines; 
new int [NUM_MACHINES] ; 


display - = new displayMatrices (NUM_TASKS, NUM_MACHINES) ; 


for(int i = 0; i < NUM_MACHINES; i++) 
{ 

MACHINE_STOP_TIME[i] = 0; 
} 


}//END CONSTRUCTOR 


[LL RRR RR KRKRKRKKKREKKKKK KKK KKK KKK KKK KEE 


//PERFORM GREEDY SEARCH 
[LL RRR RR RRKRKRKKKEKK KKK KKK KKK KKKKKEEE 
public int[][] search(int[][{] ETC, boolean{] TASK_STATUS, boolean{] 
MACHINE STATUS, 
String{][{]  TASK_ORDER, int[{] CHECKSUM, int[][] 
DEPENDENCIES, 
int[][] SCHEDULE_SUMMARY) 
{ 
int minETC; 
int MIN_ETC = 10000; 
int scheduledTask 
int scheduledMachine 
int tasks _completed = 0; 


Oo © 


= 


System.out.printin("Commence scheduling..... Ke 
System.out.printin(); 


//ENSURE ALL TASKS ARE SCHEDULED 
for(int schedule = 0; schedule < NUM_TASKS; schedule++) 


{ 
MinETC = MIN_ETC; 
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//FIND THE MINIMUM IN THE ETC ARRAY 
for(int task = 0; task < NUM_TASKS; task++) 


{ 


//ENSURE TASK HAS NOT BEEN SCHEDULED “ 
//ENSURE TASK HAS NO IMMEDIATE DEPENDENCY 
if (TASK_STATUS [task] == false && 


CHECKSUM[task] == 0) 
{ 
for(int machine = 0; machine < NUM_MACHINES; machine++) 


{ 
if { (ETC [task] [machine] + MACHINE _STOP_TIME [machine] ) < 


MInETC) 


{ 
minETC = ETC[(task] [machine] ; 


] 
scheduledTask = task; 
scheduledMachine = machine; 

} 
} 
} 
} 


//update number of tasks completed 
tasks_completed++; 


//inform user of scheduling status 


System.out.printin("Scheduling TASK #" + 

Integer.toString(scheduledTask) + 
: on MACHINE # " + 
Integer.toString(scheduledMachine) + | 
t at time i + 


Integer.toString (MACHINE_STOP_TIME[scheduledMachine] )); 


SCHEDULE_SUMMARY [schedule] [0] scheduledtTask; 

SCHEDULE_SUMMARY [schedule] [1] scheduledMachine; 

SCHEDULE_SUMMARY [schedule] [2] = 
MACHINE STOP_TIME[scheduledMachine] ; 


SCHEDULE_SUMMARY [schedule] [3] = 
MACHINE _STOP_TIME[scheduledMachine] + 


ETC [scheduledTask] [scheduledMachine] ; 


MACHINE _STOP_TIME[scheduledMachine] = 
MACHINE _STOP_TIME[scheduledMachine] + 


ETC [scheduledTask] [scheduledMachine] ; 


//update matrices 

TASK ORDER [scheduledTask] [scheduledMachine] = 
Integer.toString(schedule) ; 

TASK STATUS [scheduledTask] 

MACHINE STATUS [scheduledMachine] 


Crue; 
true; 


if(tasks_completed < NUM_TASKS) 


CHECKSUM = reviseCHECKSUM (CHECKSUM, DEPENDENCIES, scheduledTask) ; 
} ‘ 


else 


{ 
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System.out.printlin("\nScheduling complete 


} 
} 
return SCHEDULE_SUMMARY; 
}//END METHOD "search" 


[LL RRRRRRKRKEKRKEKKKEEKKEKKEKKEKKKKKKKK KS 


//REVISE THE CHECKSUM MATRIX 


[LL RRKRRREKRKKKKKKKRKEKKEKEKKEKEKREKKKKKKE 


public int[{] reviseCHECKSUM(int[] 
scheduledTask) 
{ 
if( DISPLAY MATRIX == true ) 
{ 


CHECKSUM, 


int] [] 


System.out.println ("eee RRR KEK KKK KKK KK KEK KE KM) « 
System.out.println("*** REVISE CHECKSUM() ***"); 
System. out .printIln ("eR RRR K RK KR KKK EK KKEKEM) - 


System.out.printin(); 
} 


eo # @ @ @ 


for(int second_task = 0; second_task < NUM_TASKS; 


{ 


if( DEPENDENCIES [scheduledTask] [second_task] = 


CHECKSUM [second_task]--; 
} 


if( DISPLAY_MATRIX == true ) 


display.display1D_Integer_Matrix (CHECKSUM) ; 


return CHECKSUM; 
}//END METHOD "reviseCHECKSUM" 


}// END CLASS "search_Greedy" 
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DEPENDENCIES, 


second_task++) 


int 


//creates and writes to a file using printstream object 


import java.io.*; 
import java.awt.*; : 


public class generateResultFiles 
{ 

public PrintStream ps; 

IOS os = new I0S(); 

public int NUM_TASKS; 

public int NUM_MACHINES; 


public generateResultFiles(int numtasks,int nummachines, String 
filename) 
{ 


NUM_TASKS = numtasks; 
NUM_MACHINES = nummachines; 


String file = new String("greedysummary_" + filename) ; 


Cry 4 
ps = new PrintStream(new FileOutputStream(file) ); 


} 
catch( IOException e ) { 
System.err.printin("File not opened properly\n" + 
e.toString()); 
System.exit (1); 
} 
} 


public void addSummaryHeaders () 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth(14) ; 
os.print (ps, "JOB#") ; 


os.setJustify(IOS.LEFT) ; 
os.setWidth(14); 
os.print(ps, "MACHINE#") ; 


os.setJustify(IOS.LEFT) ; 
os.setwidth(14) ; 
os.print(ps, "START TIME") ; 


os.setJustify(IOS.LEFT) ; 
os.setWidth(14) ; 
os.print(ps, "FINISH TIME") ; 


ps.printin(); 
} 


public void addScheduleSummary( int schedule_summary[][1, 
long start, 
long finish ) 


{ 
for(int task = 0; task < NUM_TASKS; task++) 


{ 


for(int displaycolumns = 0; displaycolumns < 4; displaycolumns++) 


{ 
os.setJustify (IOS.LEFT) ; 
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ee 


os.setWidth(14); 
os.print(ps,schedule_summary [task] [displaycolumns]) ; 


} 
ps.println(); 
} 
pS.printint(): 
ps.printin("Algorithm Compute Time: m+ 
’ Long.toString(finish - start) + 
" milliseconds"); 
ps.printin("Meta-task Completion Time: " + 


Integer.toString (getMetaTaskFinishTime (schedule_summary) )); 
} 


public int getMetaTaskFinishTime( int[{][] schedule_summary ) 
{ 


int meta_task_finish_time = 0; 


for(int i = 0; i < NUM_TASKS; i++) 

{ 

1f(schedule_summary[i] [3] > meta_task_finish_time) 
meta_task_finish_time = schedule_summary[i] [3]; 

} é 
return meta_task_finish_time; 

}//END METHOD "getMetaTaskFinishTime" 


public void createSummaryFile( int schedule_summary[][], 
long start, 
long finish ) 


{ 
System.out.println("Generating Result File...."); 
System.out.printlin(); 
addSummaryHeaders(); 
addScheduleSummary (schedule_summary, start, finish) ; 
} 
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public class displayMatrices implements setupConstants { 


public int NUM_TASKS; 
public int NUM_MACHINES; 


//class IOS allows text formatting similar to C++ 
IOS os = new I0S(); 


[ [RRR EERE KEEERKREEKKEKEKEKEEE KEKE KE 


//ARGUMENTED CONSTRUCTOR 


Lf FREER EERERERERARARE EERE RE RR SS 
public displayMatrices(int numtasks, int nummachines) 
{ 
NUM_TASKS 
NUM_MACHINES 
}//END CONSTRUCTOR 


numtasks; 
nummachines; 


[LL RRRRRREEKREKRREKEEEKRRKEKRERERKEEKRERERG 


//DISPLAY 1D INTEGER ARRAY - TASKS 
ff REFERS AREER RRR ER ERA ERR 
public void display1D_Integer_Matrix( int matrix[] ) 
{ 
for(int task = 0; task < NUM_TASKS; task++) 
{ 
os.setdustify(IOS.LEFT) ; 
os.setWidth (6); 
os print (System.out, "TASK#") ; 
os.setWidth (4); 
os.print (System.out, task) ; 
os.setWidth(13); 
os.printin(System.out,matrix[task]); 
}// END FOR 
System.out.println(); 
}//END METHOD "displayMatrix" 


[LR RRERRRKEEKRERERKEEKREKERKKERKEEKEKKE EK 


//DISPLAY 1D INTEGER ARRAY - MACHINES 
J [BORO IORI IO IOI RO KI IO i kk 


public void displaylD_Integer_Machine_Matrix( int matrix[] ) 


{ 
displayMachineHeaders (); 


System.out.print ("ETC Py 
for(int machine = 0; machine < NUM_MACHINES; machine++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth(13); 
os.print (System.out,matrix[machine] ); 
}// END FOR 


System.out.printin(); 
System.out.printlin(); 


}//END METHOD "display1D_Integer_Machine_ Matrix" 


[LL PER RRERREEREEREREE ERE RR ERERR EKER F ‘ 


//DISPLAY 1D BOOLEAN ARRAY 


[ [BERR EE KEKE KRREKREREREKRKE 
public void display1D_BOOL_Matrix( boolean matrix[] ) 
{ 
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f 


for(int task = 0; task < NUM_TASKS; task++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setwidth(6); 
os.print (System.out, "TASK#"); 
os.setWidth(4) ; 
os.print (System.out, task) ; 
: System.out.printlin(matrix[task] ); 
}// END FOR 
System.out.println(); 
}//END METHOD "display1D_BOOL_Matrix" 


[ [RRR RKRREKRKKKRKKRKEKKEKKKEKKKK KR KKK KK KKK 


//DISPLAY 1D STRING ARRAY 
LL RE RRRRKKRKEKRKREKEKKKKKK KKK KK KK KKK KK 
public void display1D_String Matrix( String matrix[] 
{ 
for(int task = 0; task < NUM_TASKS; task++) 
{ 
os.setdJustify(IOS.LEFT) ; 
os.setWidth(6) ; 
os.print (System.out, "TASK#") ; 
os.setWidth(4); 
os.print (System.out, task) ; 
System.out.printin(matrix[task]); 
}// END FOR 
System.out.println(); 
}//END METHOD "display1D_String_Matrix" 


[ [RRR RRRKEKEKEKKKKKKKKRKKKKKK KKK KKK 


//DISPLAY 1D FLOAT ARRAY 
f [ERP RELEASE RAR A ER ERE AERA EE ER ERE 
public void display1lD_Float_Matrix( double matrix[] ) 
{ 
for(int task = 0; task < NUM_TASKS; task++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth(6); 
os.print (System.out, "TASK#") ; 
os.setwidth(4); 
os.print (System.out, task) ; 
System.out.println(matrix[task]) ; 
}// END FOR 
System.out.printin(); 
}//END METHOD "display1D_String_Matrix" 


[LL RRR RRR KKK KKERKKK KKK KKK KKK KKK 


//DISPLAY 2D INTEGER ARRAY 
[LL RRRRRKRRERKREKEKRKEKKKKKKKKK KK KKK KKK 
public void display2D_Integer_Matrix( int matrix[][] 


{ 
displayMachineHeaders(); 


for(int task = 0; task < NUM_TASKS; task++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setwWwidth(6) ; 
os.print (System.out, "TASK#") ; 
os.setWidth(4) ; 
os.print (System. out, task); 
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) 


) 


for(int machine = 0; machine < NUM_MACHINES; machine+t+) 
{ 
os.setdJustify(IOS.LEFT) ; | ‘ 
os.setWidth(13); 
os.print (System.out,matrix[task] [machine] ) ; 
}// END FOR 
System.out.printin() ; : 
}// END FOR 
System.out.println(); 
}//END METHOD "“display2D_Integer_Matrix" 


[ [RRR KRRKREKREKRKEKREEKEKKKKEKKEKEKKEKKE 


//DISPLAY SCHEDULING ARRAY 


[LL RRR RRRRREREKEEKKKERERKEEKEKEKEKEKEKKKEEKS 
public void displaySchedulingSummary( int schedule_summary[][] ) 


{ 
/f displayMachineHeaders(); 


for(int task = 0; task < NUM_TASKS; task++) 
{ 


for(int displaycolumns = 0; displaycolumns < 4; displaycolumns++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth(13); 
os.print (System.out,schedule_summary [task] [displaycolumns]) ; 
}// END FOR 
System.out.printlin(); 
}// END FOR 
System.out.println(); 
}//END METHOD "displaySchedulingSummary" 


LL RRR RRR RRR KRKEKRKEKKEKKEKKKEKRKEKKKERS 


//DISPLAY 2D DEPENDENCY ARRAY 


LL RERRRRKERKEKRKKKEKEKKKRKKRKKKKKKKKEKER 
public void display2D_Dependency_Matrix( int matrix[][] ) 


{ 
displayTaskHeaders(); 


for(int taskl = 0; taskil < NUM_TASKS; taskl++) 
{ 

os.setJustify(IOS.LEFT) ; 

os.setWidth (6) ; 

os.print (System.out, "TASK#") ; 

os.setWidth (4); | 

os.print (System.out,taskl) ; 


for(int task2 = 0; task2 < NUM_TASKS; task2++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth(13) ; 
os.print (System.out,matrix[task1] [task2]); 
}// END FOR ; 
System.out.println(); 
}// END FOR 
System.out.printlin() ; 
}//END METHOD "display2D_Dependency_Matrix" 


[LL REREREREKREREREREKRERKREKRERKEKEREEEK 


//DISPLAY 2D STRING ARRAY 
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[ [RRR ERRKKKKRKKEKEKEKEEKKKKKKKKKKEKKKKS 
public void display2D_String_Matrix( String matrix[][] ) 


{ 
displayMachineHeaders(); 


for(int task = 0; task < NUM_TASKS; task++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth (6); 
os.print (System.out, "TASK#") ; 
os.setWidth(4); 
os.print (System.out,task) ; 


for(int machine = 0; machine < NUM_MACHINES; machine++) 
{ 
os.setdJustify(IOS.LEFT) ; 
os.setWidth(13); 
os.print (System.out,matrix[task] [machine]); 
}// END FOR 
System.out.printin(); 
}// END FOR 
System.out.println(); 
}//END METHOD "display2D_Integer_Matrix" 


[ [RRR KKRKERKEKKEREKKEKKEKRKKKEKEKER 


//DISPLAY 2D DOUBLE ARRAY 
[ [RRR KEEKKKEKEKKKKKEKKKKEKKK 
public void display2D_Double_Matrix( double matrix[][] ) 
{ 
for(int task = 0; task < NUM_TASKS; task++) 
{ 
os.setJustify(IOS.LEFT) ; 
os.setWidth(6); 
os.print (System.out, "TASK#") ; 
os.setWidth(4); 
os.print (System.out,task) ; 


for(int machine = 0; machine < NUM_MACHINES; machine++) 
{ 
os.setdustify(IOS.LEFT) ; 
os.setWidth(13); 
os.print (System.out,matrix[task] [machine] ) ; 
}// END FOR 
System.out.printin(); 
}// END FOR 
System.out.printlin(); 
}//END METHOD "displayMatrix" 


[LR RRRREKRKRKKEKKEKKRKKREKRKKEREKKKEKEKKKKE KS 


//DISPLAY MACHINE HEADERS 
[ [RRR RERKRREKREKERKRKEEKEKEKKKKKKEKEKEE 
public void displayMachineHeaders () 
{ 
System.out.print (" ae 
for(int machine = 0; machine < NUM_MACHINES; machine++) 
{ 
os.setdJustify(IOS.LEFT) ; 
os.setWidth(9); 
os.print (System.out, "MACHINE #"); 
os.setWidth(4); 
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os.print (System.out,machine) ; 

}// END FOR 

System.out.printin(); 

}// END METHOD "displayColumnHeaders" 


[ LER KKREKRKEKKEKKEKEKRKKKKKKKK KKK 


//DISPLAY TASK HEADERS 


[LER KRKKKEKKREEKKEKKEKEKKKEKK KKK 


public void displayTaskHeaders () 


{ 


System.out.print(" 
for(int task = 0; 


{ 


aa ie 


task < NUM_TASKS; task++) 


os.setdJustify(IOS.LEFT) ; 
os.setWidth (9) ; 

os.print (System.out, "TASK# "); 
os.setwidth (4); 

os.print (System.out,task); 

}// END FOR 

System.out.println(); 

}// END METHOD "displayTaskHeaders" 


[LER REKERKEKKEKEEEKEKEEKKKRKEEEKKEKKKK 


//DISPLAY INITIAL SET UP 


[LL BRRRRRKRRKEKRKRERKKERRKEREREKRKEREKEKEKERE 


public void displayInitialization(int[j[] ETC, 
DEPENDENCIES, 
int[] CHECKSUM) 
{ 
if( DISPLAY_MATRIX == true ) 
{ 
SYVSteEM- OU Dr Ine (Stee eee eee eRe eee | : 
System.out.printin("*** INITIALIZING ETC ***"); 
Svstem, Out. DELnNc ln CUTE ts ee Re ee eee R RSS ee Ree : 


} 


System. 


out.println(); 


display2D_Integer_Matrix( ETC ); 


System. 
System. 
System. 
System. 


Cut printin tl’ *** DEADLINES 


Out .printin(); 


display1lD_Integer_Matrix( DEADLINE ); 


System. 
System. 
System. 
System. 


out.println("*** DEPENDENCIES 


out .printlin(): 


display2D_Dependency_Matrix( DEPENDENCIES ); 


SVYSteM.ouE<DrinNeini(" =e Asse ee See Ae Ee ee) ; 
kKeEN) ; 


System 
System 
System 


.OUE .<DrIntan(**** CHECK SUM 


SOUT.  PLintIn (TEAK KHER EER KEKE ARE EM) 5 


.out.println(); 


display1lD_Integer_Matrix( CHECKSUM ); 


}//END METHOD "displayiInitialization" 


}//END CLASS "displayMatrices" 
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SUL ADEN LCT tea ee eR EE ee ee) : 
mee : 


SUL DEINE li ses st FAs eee eA ee SARS) : 


OUC .PLint ln (UKE K HRA HHRKEKKKKKEKEEEN) 5 
weEM Ys 


Cue Seine Init ee Akt ee eee eee eee ee eee 


DEADLINE, 


int[][] 








) a 


* Frederick F. Chew, The Java/C++ Cross Reference Handbook (Book/CD- 
ROM) 
* Published By HP Press/Prentice-Hall 
* Copyright (C) 1997 Hewlett-Packard Company 
* All Rights Reserved. ISBN 0-13-848318-3 
* 
* Permission to use, copy, modify, and distribute this 
* software and its documentation for NON-COMMERCIAL purposes 
* and without fee is hereby granted provided that this 
* copyright notice appears in all copies. 
* 
* THE AUTHOR AND PUBLISHER MAKE NO REPRESENTATIONS OR 
* WARRANTIES ABOUT THE SUITABILITY OF THE SOFTWARE, EITHER 
* EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE 
* IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A 
* PARTICULAR PURPOSE, OR NON-INFRINGEMENT. THE AUTHOR 
* AND PUBLISHER SHALL NOT BE LIABLE FOR ANY DAMAGES SUFFERED 
* BY LICENSEE AS A RESULT OF USING, MODIFYING OR DISTRIBUTING 
* THIS SOFTWARE OR ITS DERIVATIVES. 
es 
[RX 
* @version 1.00 01 May 1997 
* @author Frederick F. Chew 
Be 4 
[** 
* Note This program was built from the JDK version 1.02 
= and has not been customized to work on version 1.1 
* 
- The precision of floating point numbers has undergone 
: change in version 1.1 
al 
import java.awt.*; 
import java.io.*; 
public class I0S { 
public I0S() { 
radix = DEC; 
width = 0; 
El iechar. = ms «3 
justify_right = RIGHT; 
precision = 6; 
upper_case = false; 
show_sign = false; 
internal pad = false; 
fixed_format = FIXED; 
} 
public IO0S(int base, int size, char fill, boolean just, int 


Significant, boolean cap, 
boolean sign, boolean pad, boolean format) { 


if ((base < 0) || (base > 2)) 
base = DEC; 
if (size < 0) 


Size = 0Q; 
radix = base; 
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width = size; 

£7 LGnar Ss: Paid bs: 

justify_right = just; 

precision. = -significant; ; 
upper_case = cap; 

show_sign = sign; 

internal_pad = pad; 

fixed_format format; 


} 


// public mutator methods for IOS attributes 


public int setRadix(int base) { 
int org_radix = radix; 
if ((base < 0) || (base > 2)) 
base = DEC; 
radix = base; 
return org_radix; 


} 


public int setWidth(int size) { 
int org_width = width; 
if (size < 0) 
size = Q; 
width = size; 
return org_width; 


} 


public char setFill(char fill) { 
char org_fillchar = fillchar; 
fillchar = £111; 
return org_fillchar; 


} 


public boolean setJustify(boolean just) { 
boolean org_just = justify_right; 
justify_right = just; 
return org_just; 


} 


public int setPrecision(int prec) { 
int org _precision = precision; 
precision = prec; 
return org _precision; 


} 


public boolean setUpperCase(boolean cap) { 
boolean org_upper_case = upper_case; 
upper_case = cap; 
return org_upper_case; 


} 


public boolean setShowSign(boolean sign) { 
boolean org _show_sign = show_sign; 
show_sign = sign; 
return org_show_sign; 


} 


public boolean setInternalPad(boolean pad) { 
boolean org_internal_pad = pad; 
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fo 


internal_pad = pad; 
return org_internal_pad; 


} 


public boolean setNumericalFormat (boolean format) { 
boolean org_fixed_format = fixed_format; 
fixed_format = format; 
return org_fixed_format; 


} 


void showState(PrintStream ps) { 
String: radix ser =--'; 
String width_str = ""; 
String just_str = ""; 
String fillchar_str = ""; 
String format_str = ""; 





Switch (radix) { 
case DEC: radix_str += "Decimal"; 


break; 

case OCT: radix str += "Octal"; 
break; 

case HEX: radix_str += "Hexadecimal"; 
break; 


} 


if (width == 0) 
width_str += "Default"; 
else 

width str += width; 


if (justify_right) 
just_str += "RIGHT"; 
else 
just_str += "LEFT"; 


LG tridlichar == ~~  *) 
fillchar_str += "Space"; 
else 
fillenar str +=] fillenar: 


if (£fixed_format) 
format_str += "Fixed"; 
else 
format_str += "Scientific"; 


ps.println(); 


ps.println("Radix = "+radix_str) ; 
ps.printlin("Width = "+width_str) ; 
ps.println("Fill Character = "+fillchar_str) ; 
ps.printlin("Justification = "+just_str); 
ps.println("Precision = "+precision) ; 
ps.println("Upper case = "+upper_case); 
ps.printlin("Show sign = "+show_sign) ; 
ps.printlin("Internal pad = "+internal_pad) ; 

, ps.printin("Numerical format = "+format_str); 
ps.println(); 

} 


// public methods for print() and printin() 
51 





public void println(PrintStream ps, char aChar) { 
String newStr = ""; 


newStr += aChar; $ 
print(ps, newStr) ; 
ps.println(); 


} 


public void printlin(PrintStream ps, char [] aCharArray) { 
String newStr = ""; 
newStr += aCharArray; 
print(ps, newStr) ; 
ps.printin(); 
} 


public void printin(PrintStream ps, int anInt) { 
println(ps, (long)anInt) ; 
} 


public void println(PrintStream ps, long aLong) { 
print(ps, aLong); 
PS.printin(); 

} 


public void printin(PrintStream ps, float aFlt) { 
println(ps, (double) aFlt); 
} 


public void printlin(PrintStream ps, double aDbl) { 
print(ps, aDbl); 
ps.println(); 

} 


public void printin(PrintStream ps, String aStr) { 
print(ps, aStr); 
ps.printin(); 

} 


public void print (PrintStream ps, char aChar) { 
String newStr = ""; 
newStr += aChar; 
print(ps, newStr) ; 


} 


public void print(PrintStream ps, char [] aCharArray) { 
String newStr = ""; 
newStr += aCharArray; 
print(ps, newStr); 


} 


public void print (PrintStream ps, int anInt) { 
print(ps, (long)anInt) ; 
} 


public void print (PrintStream ps, long aLong) { 
String newStr = ""; 
1f ((show_sign) && (along > 0)) 
newStr += "+"; 
switch (radix) { 
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case DEC: newStr += abLong; 
break; 


case OCT: newStr += "0" + newBase(aLong, 3, 7, "01234567"); 


break; 
case HEX: if (upper_case) 
newStr += ~O%" + newBase (aLong, 4, 
"0123456789ABCDEF") ; 
else 
newStr += Ose" + newBase (aLong, 4, 
"0123456789abcdef") ; 
break; 
} 
print(ps, newStr); 
} 


public void print(PrintStream ps, float aFlt) { 
print(ps, (double)aFlt) ; 
} 


public void print(PrintStream ps, double aDbl) { 
String newStr = ""; 
String tmpStr = ""; 
if (fixed_format) 
tmpStr += prepareFixedFormat (aDbl1) ; 
else 
tmpStr += prepareScientificFormat(aDbl) ; 
newStr += padString(tmpStr, aDbl); 
print(ps, newStr); 


} 


public void print(PrintStream ps, String aStr) { 
ps.print (prepareField(aStr)); 
ps.flush(); 

} 


// private helper methods 


private void drawRuler(PrintStream ps) { 
String interval = new String("1234567890") ; 
String ruler = ""; 
for (int index = 0; index < 6; index++) 

ruler += interval; 

pS.printin(); 
ps.printin(ruler) ; 
ps.printin(); 

} 


sy 


d Bo 


private String newBase(long value, int no_of bits, int one_bits, 


String digits) { 

String result = ""; 

if (value == Q) { 
result += "0"; 

} 

else { 
while (value != 0) { 
result = digits.charAt((int) (value & one_bits)) + result; 
value = value >>> no_of_bits; 
} 

} 


return result; 
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} 


private String prepareField(String aStr) { 
if. (Wreath == ..0): 4 
return aStr; 


} 


if (width <= aStr.length()) { 
return aStr; 


} 


un. 
é 


String newStr = 
int delta = width - aStr.length(); 
if (justify_right) { // RIGHT 
for (int index = 0; index < delta; 
newStr += fillchar; 
newStr += aStr; 
} 
else { 
StringBuffer sbuffer = 
char [] cbuffer new char[delta]; 
sbuffer.append(cbuffer) ; 
for (int index = 0; index < delta; 


// LEFT 


sbuffer.setCharAt (aStr.length()+index, 


newStr += new String(sbuffer) ; 


} 


return newStr; 


} 


private String eNotation(char sign, 
String expStr = 
int leading_zeros = 3 - 
String result = ""; 
result += (upper_case) ?"E":"e"; 
result += sign; 
for (int index = 

result += "Q"; 

result += expStr; 
return result; 


} 


0; 


private String moveDecimalLeft(long aLong) 
String result = ""; 
String Shes: Ve. 


str += aLong; 
result += str.substring(0, 1) + 
str.length()); 


return result; 


} 


private String moveDecimalRight (String str) 
char ch: 
long pwr = 0; 
int index2 = 0 
String result 
String newStr 


indexl < str.length(); 
Oe 3 


(int indexl = 0; 
((ch = str.charAt (index1) ) 


ror 
if 
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index < leading zeros; 


index++) 


new StringBuffer(aStr) ; 


index++) 


fac lchar):; 


long exponent) {f{ 
new Long(exponent) .toString(); 
(expStr.length()),; 


index++) 


{ 


+ str.substring(1, 


{ 


indexl++) { 








PWIr++t ; 
else 
newStr += ch; 
: } 
PWIr++; 
» if ((newStr.substring(1, newStr.length())).length() != 0) 
result += newStr.substring(0, 1) + "." + #£=newStr.substring(1, 
newStr.length()) + eNotation(’-’, pwr); 
else | 
result += newStr.substring(0, 1) + "." + "0" + eNotation(’-’, pwr); 


return result; 
} 


private boolean alreadyScientific(String str) { 

boolean flag = false; 
for(int index=0; index < str.length(); index++) { 

if ( (str.charAt (index) == ‘e’) || (str.charAt(index) == 'E’) ) 
flag = true; 

} 
return flag; 
} 





private String improveDefault(String str) { 
boolean decimal_pt = false; 
int e_index = 0; : 
for (int index = 0; index < str.length(); index++) { 





Lf (StriecnearaAt( index). == 7.7) 
decimal_pt = true; 
else 
if ( (str.charAt (index) == ’e’) || (str.charAt(index) == ’E’) ) 


e index = index; 


} 
String result = ""; 


if (decimal_pt) { 





result += str.substring(0, e_index) + (upper_case?"E":"e") + 
str.substring(e_index+i, str.length()); 
else { 
result += str.substring(0, e_index) + ".0" + (upper_case?"E":"e") 
+ str.substring(e_index+1, str.length()); 





} 


return result; 


} 


private String padString(String str, double num) { 
int lgth = width - str.length(); 


if (((show_sign) || (num < 0.0)) && (lgth > 0) && (internal_pad)) { 
String result = ""; | 
result += str.substring(0, 1); 
for (int index = 0; index < lgth; index++) result += fillchar; 
result += str.substring(l, str.length()); 
return result; 
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return str; 


} 


private String prepareScientificFormat (double aDbl) { 
String result = ""; 


// The double may already be in scientific format 


String dblStr = ""; 
ab1lStr += aDbl; 


if (alreadyScientific(dblStr)) { 
String newStr = ""; 
if ((show_sign) && (aDbl > 0)) 
newStr += "+"; 
newStr += improveDefault (db1lStr) ; 
return newStr; 


} 


// Include the sign if necessary 


if (aDbl < 0.0) { 
result += "-"; 
aDbl = -1 * aDbl; 
J 


else if (show_sign) 
result += "+"; 


// Find the whole part and the fractional part 


long whole = (long) aDbl; 
double fraction = aDbl - whole; 


String wholeStr = ""; 
wholeStr += whole; 


if (wholeStr.length() >= precision) { 
if (fraction >= 0.5) wholet++; 
long pwr = wholeStr.length() - 1; 
result += moveDecimalLeft (whole) + eNotation(’+’, pwr); 


} 


else { 
int frac_prec = precision - wholeStr.length(); 
if (fraction == 0.0) { 
long pwr = wholeStr.length() - 1; 
result += moveDecimalLeft (whole) + ((whole==0) ?"0":"") + 


eNotation(’+’, pwr); 
} 


else 
result += determineFraction(frac_prec, fraction, whole); 


} 


return result; 


} 


private String determineFraction(int prec, Gouble aDbl, long 
org_whole) { 


String result = ""; 
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String abblstr-= ""; 

aDblStr += aDbl; 

int lgth = aDblStr.length({)-2; 
Strang frac.str: =" ™"; 


if (prec < lgth) { 
double newDbl = aDbl; 
for (int index = 0; index < prec; index++) 
newDbl *= 10.0; 
long sub_whole = (long)newDbl; 
double sub_fraction = newDbl - sub_whole; 
1£ (sub_fraction >= 0.5) sub_whole++; 
frac_str += sub_whole; 
} 
else 
frac_str += aDblStr.substring(2, aDblStr.length()); 


String orgWholeStr = ""; 
orgwholeStr += org_whole; 


if (org _whole > 0) { 
long pwr = orgWholeStr.length() - 1; 
result += moveDecimalLeft(org._whole) + frac_str + eNotation(’+’, 
pwr); 
} 
else { 
result += moveDecimalRight (frac_str) ; 


i 


return result; 


} 


private String prepareFixedFormat (double aDbl) { 
String result = "2 


// The double may already be in scientific format 


String dblistr: =-"" 
dblStr += aDbl; 


1£ (alreadyScientific(dblStr)) { 
String newStr = ""; 

if ((show_sign) && (aDbl > 0)) 
newStr += "+"; 

newStr += improveDefault(dblStr) ; 

return newStr; 


} 
// Include the sign if necessary 


if (aDbl < 0.0) { 
result += "-"; 
aDbl = -1 * aDbl; 

} 

else if (show_sign) 
result += "+"; 


// Find the whole part and the fractional part 


long whole = (long) aDbl; 
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double fraction = aDbl - whole; 


String wholeStr = ""; 
wholeStr += whole; 


L£ (wholeStr.length() >= precision) { 
if (fraction >= 0.5) whole++; 
result += whole; 


} 


else { 
int frac_prec = precision - wholeStr.length(); 
if (fraction == 0.0) 
result += whole; 
else 
result += whole + "." + findFraction(frac_prec, fraction) ; 
} 


return result; 


} 


private String findFraction(int prec, double aDbl) { 
String result = ""; 

String aDblStr = ""; 

aDbiStr += aDbl; 

int lgth = aDblStr.length()-2; 


if (prec < lgth) { 
double newDbl = aDbl; 
for (int index = 0; index < prec; index++) 
newDbi *= 10.0; 
long sub_whole = (long)newDbl; 
double sub_fraction = newDbl - sub_whole; 
if (sub_fraction >= 0.5) sub_whole++; 
result += sub_whole; 
} 
else 
result += aDblStr.substring(2, aDblStr.length()); 


return result; 


// public final fields 


public static final int DEC = 
public static final int OCT = 1; 
public static final int HEX = 2; 


public static final boolean RIGHT = true; 

public static final boolean LEFT = false; 
public static final boolean FIXED = true; 

public static final boolean SCIENTIFIC = false; 


// private variable fields 


private int precision; 
private int radix; 

private int width; 

private char fillchar; 

private boolean justify_right; 
private boolean upper_case; 
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private boolean show_sign; 
private boolean internal_pad; 
private boolean fixed_format; 


‘ public static void main(String args[]) throws IOException { 


double [] dbli = { 

D205: YOG-1s 02013 00012 20-0001: 02 00001, 

0.000001, 0.0000001,  0.00000001, 0.000000001, 0.0000000001, 
| 0.00000000001 


ie 


double [] dbl2 = { 
2 kee ie Oi Ze Ui 254 0.1234; Or i2545. 0.123456, 0.1234567, 
0.12345678, 0.123456789, 
0.1234567891, 0.12345678912, 0.123456789123, 0.1234567891234, 
0.12345678912345 
iG 


double [] dbl3 = { 





0.123456, Lat Z2s256; 124123456; £23 .123456; 1234.123456, 
12345 .123456,. 123456.:123456, 

12345 .6123456, 1234.56123456;, 123 .456123456, 12.3456123456, 
1.23456123456, 

1234561.23456, 12345612 .3456, 123456123 .456, 1234561234.56, 
12345612345.6, 1234567.123456,. -12345678.123456; 
123456789 .123456, 1234567891 .123456, 12345678912.123456, 
123456789123 .123456, 1234567891234 .123456 

| Ps 


double [] dbl4 = { 
cee Ieee T2412; 123.123, 1234.1234, 12345. 1l2540% 123456.123456, 


1234567 .1234567, 12345678 .12345678, 123456789 .123456789, 
1234567891 .1234567891, 12345678912 .12345678912, 
123456789123 .123456789123, 

1234567891234 .1234567891234, 12345678912345 .12345678912345, 
123456789123456 .123456789123456, 1234567891234567 .1234567891234567, 


12345678912345678 .12345678912345678, 
123456789123456789 .123456789123456789 
ie 


// TestoOl 
PrintStream prsl = new PrintStream (new 


FileOutputStream("OS1A.TXT")); 
IOS oslA = new IOS(); 





String bird = "penguin"; 
String fish = "rock cod"; 
String tree = "maple"; 


‘ oslA.setFill(’~’); 
oslA.setWidth(12); 
osiA.printin(prsl, bird); 
oslA.printin(prsl, fish); 
oslA.setWidth (0) ; 
osiA.println(prsl, tree); 


prsi.close(); 





// Test02 


PrintStream prs2 = new PrintStream (new 


FileOutputStream("OS2A.TXT") ); 
TOS os2A = new IOS(); 


double numl = 17.7875678; 
os2A.setPrecision (8); 
os2A.println(prs2, numl); 


prs2.close(); 


// Test03 


PrintStream prs3 = new PrintStream(new 


FileOutputStream("OS3A.TXT")); 
IOS os3A = new IOS(); 


int num2 = 160; 


os3A.setRadix(I0OS.OcT) ; 
os3A.print(prs3, "Octal = "); 
os3A.println(prs3, num2); 
os3A.setRadix (IOS.HEX) ; 
os3A.print(prs3, "Hexadecimal = "); 
os3A.printin(prs3, num2); 
os3A.setUpperCase (true) ; 
os3A.print(prs3, “Hexadecimal = "); 
os3A.println(prs3, num2) ; 
os3A.setRadix(IOS.DEC) ; 
os3A.print(prs3, "Decimal = "); 
os3A.printin(prs3, num2); 


prs3.close(); 


// Testod4 


PrintStream prs4 = new PrintStream(new 


FileOutputStream("OS4A.TXT")); 
IOS os4A = new IOS(); 


“bod. cat"; 
757905; 
2345070.00915756; 


String feline 
long inum 
double dnum 


os4A.setJustify(IOS.LEFT) ; 
os4A.setFill(’#’); 
os4A.setWidth(12); 
os4A.setPrecision (14); 
os4A.printin(prs4, feline) ; 
os4A.printlin(prs4, linum); 
os4A.printin(prs4, dnum); 


os4A.setJustify (IOS.RIGHT) ; 
os4A.setFill(’*’); 
os4A.setWidth(12) ; 
os4A.setPrecision (14) ; 
os4A.println(prs4, feline); 
os4A.println(prs4, lnum); 
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os4A.println(prs4, dnum); 
prs4.close(); 
// Test05 


PrintStream prs5 = new PrintStream (new 
. FileOutputStream("OSSA.TXT") ); 
IOS osSA = new I0S(); 


ossA.setdJustify(IOS.LEFT) ; 
osSA.setPill(’*’); 
osSA.setWidth(12); 
ossSA.printlin(prs5, bird); 
os5A.printin(prs5, fish); 
ossSA.setJustify(IOS.RIGHT) ; 
osSA.println(prs5, tree); 


prs5.close(); 
// Test06 
PrintStream prs6 = new PrintStream (new 


FileOutputStream("OS6A.TXT") ); 
TOS os6A = new IO0S(); 


ny okies intl = 34950; 
int int2. = Hints 
Ficat fitl. = 17-5829; 
float flt2 =--£Lel1; 


os6A.setinternalPad (true) ; 
os6A.setFill('%’); 
os6A.setWidth(16) ; 
os6A.printin(prs6, intl); 
os6A.printin(prs6, int2); 
os6A.printin(prs6, fitl); 
os6A.printin(prs6, £1t2); 
os6A.setiInternalPad(false) ; 
os6A.printin(prs6, f1t2); 


prs6.close(); 
// Test0O7 
PrintStream prs7 = new PrintStream (new 


FileOutputStream("OS7A.TXT") ); 
IOS os7A = new IOS(); 


long longl = 875090406; 
double doubl = 395.0; 
double doub2 = 7693.8356; 
double doub3 = 3409.958; 


os7A.println(prs7, doubl) ; 
os7A.setShowSign (true) ; 
os7A.println(prs7, doubl); 
os7A.setNumericalFormat (IOS.SCIENTIFIC) ; 
os7A.println(prs7, longl); 
os7A.printin(prs7, doub2) ; 
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os7A.printin(prs7, doub3); 

os7A.setShowSign (false) ; 

os7A.setNumericalFormat (IOS.FIXED) ; 

os7A.printin(prs7, longl); , 
os7A.println(prs7, doubl); 

os7A.printin(prs7, doub2) ; 

os7A.println(prs7, doub3); 


prs7.close(); 


// Test08 


PrintStream prs8 = new PrintStream(new 


FileOutputStream("OS8A.TXT")); 
IOS os8A = new IOS(); 


os8A.setPrecision(20); 
os8A.setWidth (30); 
os8A.setJustify(IOS.RIGHT) ; 
os8A.drawRuler(prs8) ; 


os8A.printlin(prs8, "Default System.out.println() - Positive"); 
for(int index = 0; index < dbll.length; index++) 
prs8.println(dbli [index] ); 


os8A.println(prs8, "Default System.out.printin() - Negative"); 
for(int index = 0; index < dbli.length; index++) 
prs8 .printin(-dbli[index]); 


os8A.printin(prs8, "Fixed Format - Positive"); 
for(int index = 0; index < dbli.length; index++) 
os8A.printlin(prs8, dbli[index]); 


os8A.printlin(prs8, "Fixed Format - Negative"); 
for(int index = 0; index < dbil.length; index++) 
os8A.printin(prs8, ~dblifindex] ) ; 


os8A.setNumericalFormat (IO0S.SCIENTIFIC) ; 


os8A.printlin(prs8, "Scientific Format - Positive"); 
for(int index = 0; index < dblil.length; index++) 
os8A.printlin(prs8, dblif{index]); 


os8A.println(prs8, "Scientific Format - Negative"); 
for(int index = 0; index < dbli.length; index++) 
os8A.printin(prs8, -dbli[index]); 


prs8.close(); 


// Testd09 


PrintStream prs9 = new PrintStream (new 


FileOutputStream("OS9A.TXT")); 
IOS os9A = new IOS(); 


os9A.setPrecision(20) ; 
os9A.setWidth(30); 
os9A.setJustify(IOS.RIGHT) ; 
os9A.drawRuler(prs9); 
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os9A.printlin(prs9, "Default System.out.println() 
for(int index = 0; index < dbl2.length; index++) 
prs9.printiln(dbl2 [index] ); 


os9A.printin(prs9, "Default System.out.printin( 
for(int index = 0; index < dbl2.length; index++ 
prs9.printin(-dbl2 [index] ) ; 


) 


os9A.printin(prs9, "Fixed Format - Positive") ; 
for(int index = 0; index < dbl2.length; index++) 
os9A.printin(prs9, dbl2 [index] ); 


os9A.printin(prs9, "Fixed Format - Negative"); 
for(int index = 0; index < dbl2.length; index++) 
os9A.printlin(prs9, -dbl2[index]) ; 


os9A.setNumericalFormat (IOS.SCIENTIFIC); 





- Positive"); 


~ Negative"); 


os9A.printlin(prs9, "Scientific Format - Positive"); 


for(int index = 0; index < dbl2.length; index++) 
os9A.println(prs9, dbl2[index] ); 


os9A.println(prs9, "Scientific Format - Negative") ; 


for(int index = 0; index < dbl2.length; index++) 
os9A.printin(prs9, -dbl2[index]); 


ors9.close(); 
// Testl10 


PrintStream prsl0 = new 


FileOutputStream("OS10A.TXT") ); 


IOS os1l0A = new I0S(); 


osl0A.setPrecision(20); 
osl0A.setWidth(30); 
osl0A.setJustify(IOS.RIGHT) ; 
osl0A.drawRuler(prs10); 


PrintStream(new 


osl0A.println(prs10, "Default System.out.printin() - Positive") ; 


for(int index = 0; index < dbl3.length; index++) 
prsl0.println(dbl13 [index] ) ; 


osl0A.printlin(prs10, "Default System.out.println() - Negative"); 


for(int index = 0; index < dbl3.length; index++) 
prs10.println(-db1l3 [index]); 


osl0A.println(prs10, "Fixed Format - Positive") ; 
for(int index = 0; index < dbi3.length; index++) 
osl0OA.printlin(prsi0, db1l3 [index] ) ; 


osi0A.printin(prs10, "Fixed Format - Negative"); 
for(int index = 0; index < dbl3.length; index++) 
osl0A.printlin(prs10, -dbi3{findex]); 


osl0A.setNumericalFormat (IOS.SCIENTIFIC) ; 


osl0A.printlin(prsl10, "Scientific Format - Positive"); 


for(int index = 0; index < dbl3.length; index++) 
osl0A.printlin(prsl10, dbl3[index]); 
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osl0A.println(prsi0, "Scientific Format - Negative"); 
for(int index = 0; index < dbl3.length; indext++) 
osl0A.println(prs10, -db13[index]); : 


prs10.close(); 


// Testil 


PrintStream prsil = new PrintStream (new 


FileOutputStream("OS11A.TXT") ); 


} 


} 


TOS osl1lA = new IOS(); 


osllA.setPrecision(20); 
osll1A.setWidth (30); 
osliA.setJustify(IOS.RIGHT) ; 
osllA.drawRuler(prs11) ; 


osilA.println(prs11, "Default System.out.printin() - Positive"); 
for(int index = 0; index < dbl4.length; index++) 
prsil.printin(dbl4[index]); 


osilA.printlin(prsil, "Default System.out.println() - Negative"); 
for(int index = 0; index < dbl4.length; index++) 
prsil.printin(-dbl4[index]); 


osilA.printin(prsll1, "Fixed Format - Positive"); 
for(int index = 0; index < dbl4.length; index++) 
osllA.printin(prsil, dbl4[index]); 


osl1A.println(prs11, "Fixed Format - Negative"); 
for(int index = 0; index < dbl4.length; index++) 
osliA.println(prsii, -dbl14[index] ); 


osliA.setNumericalFormat (IOS.SCIENTIFIC) ; 


osllA.println(prsi11l, "Scientific Format - Positive"); 
for(int index = 0; index < dbl4.length; index++) 
osllA.println(prsl11, dbl4[index]); 


osilA.println(prs11, "Scientific Format - Negative") ; 
for(int index = 0; index < dbl4.length; index++) 
osliA.printlin(prsil, -dbl4[index] ); 


prsli.close(); 


System.out.println("(press Enter to exit)"); 
LY 4 
System.in.read(); 


} 
catch (IOException e) { 


return; 


} 
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public interface setupConstants { 


public 
public 
public 
public 
public 
public 
public 


public 
public 


int dependency_type = RANDOM_DEPENDENCY; 


finak 
final 
final 
final 
final 
final 
final 


final 
final 


Static int ETC RANGE 


static boolean DISPLAY MATRIX 
static boolean DISPLAY _TASKS 
static int FULL_DEPENDENCY 


static int NO_DEPENDENCY 
static int RANDOM DEPENDENCY 


static int GREEDY 
static int A_STAR 
static int GENETIC 


int search_type = GREEDY; 


}//END INTERFACE "setupConstants" 
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ri; astar search algorithm 


(defun set-job-list () 
(let* ((job-list nil)) : 
(dotimes (i *number-of-jobs* job-list) 
(setf job-list 
(append job-list (list i)))))) 


(defun set-machine-list () 
(let* ((machine-list nil)) 
(dotimes (i *number-of-machines* machine-list) 
(setf machine-list 
(append machine-list (list i)))))) 


(defun get-time-to-compute (job machine) 
(elt (elt *etc-matrix* job) machine) ) 


(defun astarsearch (start input-file output-file) 
(read-input-file input-file) 
(setf *algorithm-start-time* (get-universali-time) ) 
(set£ *schedule* (astarsearch2 
(list (list (evalfn start) start)) nil)) 
(setf *algorithm-stop-time* (get-universal-time) ) 
(write-output-file output-file) ) 


(defun astarsearch2 (agenda usedagenda) 
(when agenda 
(2E 

(goalreached (second (first agenda) )) 
(last (reverse (rest (first agenda) ))) 
(astarsearch2 

(merge 

eTESt 

(remove-bad-agitems 

(mapcar 
#’ (lambda 
(newstate) 
(let ((newpath 
(cons newstate (rest (first agenda))))) 
(cons (+ (evalfn newstate) (costfin newpath) ) 
newpath) )) 
(successors (second (first agenda) ))) 
(append agenda usedagenda)} ) 

(rest agenda) 

#’firstlessp) 

(cons (first agenda) usedagenda) )))) 


(defun remove-bad-agitems (agitems oldagitems) 

(sort 

(remove-if 

#’ (lambda (agitem) F 
(some #’ (lambda (oldagitem) 
(and (>= (first agitem) (first oldagitem) ) 
(equalstate (second agitem) (second oldagitem)))) 
oldagitems) ) 
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agitems) 
#‘firstlessp) ) 


(defun firstlessp (agiteml agitem2) 
(< (first agiteml) (first agitem2) )) 


(defun successors (state) 
(let* ((immediate-successors (successors2 state) ) 
(complete-successors nil)) 
(dotimes (1 (length immediate-successors) complete-successors) 
(setf complete-successors 
(append complete-successors 
(list (append (elt immediate-successors 1) 

state))))))) 


{(defun successors2 (state) 
(let* ((successors nil) 
(jobs-to-schedule (get-jobs-to-schedule state) ) 
(machine-stop-times (get-machine-stop-times state) )) 
(dotimes (m (length *machines*) successors) 
(dotimes (j (length jobs-to-schedule) nil) 
(1f (dependency-status (elt jobs-to-schedule j) 
jobs-to-schedule) 





(setf successors 
(cons 
(list (list (elt jobs-to-schedule j) 
(elt *machines* m) 
(elt machine-stop-times m) 
(+ (elt machine-stop-times m) 
(get-time-to-compute (elt sjobs-to-schedule 


successors))))))) 


(defun goalreached (statelist) 
(if (equalp (length (get-jobs-to-schedule statelist)) 0) t nil)) 


(defun equalstate (statel state2) 
(equal statel state2)) 


(defun costin (statelist) 
(let* ((cost 0) 
(machine-cost (get-individual-machine-cost (first . 
statelist)))) 
(dotimes (n (length machine-cost) cost) 
(1£ (> (elt machine-cost n) cost) 
(sett cost (elt machine-cost n)))))) 





(defun evalfn (statelist) 
(if (= (length (get-jobs-to-schedule statelist)) 0) 0 
(1£ (S= (length (get-jJobs-to-schedule statelist) ) (length 
*machines*) ) 
(/ (get-unscheduled-job-time statelist) (length *machines*) ) 
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(/ (get-unscheduled-job-time statelist) (length (get-jobs-to- 
schedule statelist)))))) 


(defun get-jobs-to-schedule (state) 
(let* ((jobs-to-schedule *jobs*) ) 
(dotimes (n (length state) jobs-to-schedule) 


(if (memberofp (first (elt state n)) jobs-to-schedule) , 
(set£ jobs-to-schedule 
(remove (first (elt state n)) jobs-to-schedule)))))) 


(defun get-average-job-time () 
(let* ((average-job-time (init-job-list) )) 
(dotimes (i (length *jobs*) average-job-time) 

(dotimes (j (length *machines*) ) 

(setf (elt average-job-time i) 
(+ (elt average-job-time i) 
(get-time-to-compute i j)))) 

(set£f (elt average-job-time 1) 

(/ (elt average-job-time i) (length *machines*)))))) 


(defun get-unscheduled-job-time (statelist) 
(let* (({time-left 0) 
(jobs-remaining (get-jobs-to-schedule statelist) )) 
(dotimes (n (length jobs-remaining) time-left) 
(setf time-left 
(+ time-left (elt (get~average-job-time) (elt jobs-remaining 
n))))))) 


(defun get-individual-machine-cost (state) 
(let* ((machine-cost (init-machine-list) )) 
(dotimes (n (length state) machine-cost) 
(let* ((machine-to-update (second (elt state n))) 
(job-start-time (third (elt state n))) 
(job-finish-time (fourth (elt state n))) 
(job-scheduled (first (eit state n)))) 
(setf (elt machine-cost machine-to-update) 
(+ (elt machine-cost machine-to-update) 
(get-time-to-compute (first (elt state n)) 
(second (elt state n))))) 
(if (> job-finish-time (elt *deadline-matrix* job-scheduled) ) 
(set£ (elt machine-cost machine-to-update) 
(+ (elt machine-cost machine-to-update) 
(- job-finish-time Celt *deadline-matrix* job- 
scheduled))))))))) 


(defun get-machine-stop-times (state) 
(let* ((machine-stop-times (init-machine-1list) )) 
(dotimes (i (length state) machine-stop-times) > 
(let* ((machine (second (elt state i)))) 
(setf (elt machine-stop-times machine) 
(fourth (elt state i))))))) : ; 


(defun init-job-list () 
(let* ((average-job-time nil)) 
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ee 


(dotimes (i (length *jobs*) average-job-time) 
(setf average-job-time 
{append ’(0) average-job-time))))) 


(defun init-machine-list () 
(let* ((machine-cost nil)) 
‘ (dotimes (i (length *machines*) machine-cost) 
(setf machine-cost 
(append ‘(0) machine-cost))))) 


(defun notmemberofp (item list) 
(let* ((memberstatus t)) 
(dotimes (n (length list) memberstatus) 
(1£ (equalp item (elt list n)) 
(setf memberstatus nil))))) 


| (defun memberofp (item list) 
(let* ((memberstatus nil)) 
(dotimes (n (length list) memberstatus) 
(1£ (equalp item (elt list n)) 
(setf memberstatus t))))) 





(defun read-input-file (filename) 
(setq in-stream (open filename :direction :input)) 
(let* ((number-of-jobs (read in-stream) ) 
(number-of-machines (read in-stream) ) ) 

(set£f *number-of-jobs* number-of-jobs) 
(sett *number-of-machines* number-of-machines) 
(setf *jobs* (set-job-list) ) 
(setf *machines* (set-machine-list) ) 
(read-etc-matrix in-stream) 
(read-dependency-matrix in-stream) 
(read-deadline-matrix in-stream) 
(close in-stream) ) ) 


(defun read-etc-matrix (in-stream) 
(let* ((individual-job-times nil) 
(etc-matrix nil) ) 
(dotimes (i *number-of-jobs*) 
(set£ individual-job-times nil) 
(dotimes (j *number-of-machines*) 
(setf£ individual-job-times 
(append individual-job-times (list (read in-stream) )))) 
(set£ etc-matrix 
(append etc-matrix (list individual-job-times) )) ) 
(setf£f *etc-matrix* etc-matrix) )) 


(defun read-dependency-matrix (in-stream) 
(let* ((individual-job-dependencies nil) 
(dependency-matrix nil) ) 
(dotimes (1 *number-of-jobs*) 
(setf individual-job-dependencies nil) 
(dotimes (j *number-of-jobs*) 
(setf individual-job-dependencies 
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{append individual-job-dependencies (list (read in- 


stream) )))) 
(setf dependency-matrix 
(append dependency-matrix (list individual-job- : 
dependencies) ))) 


(set£ *dependency-matrix* dependency-matrix) ) ) 


(defun read-deadline-matrix (in-stream) 
(let* ((deadline-matrix nil)) 
(dotimes (i *number-of-jobs*) 
(set£ deadline-matrix 
(append deadline-matrix (list (read in-stream))))) 
(setf *deadline-matrix* deadline-matrix) ) ) 


(defun write-output-file (filename) 
(setq out-stream (open filename :direction :output) ) 
(print ‘LISP-A*-Search-Summary out-stream) 
(dotimes (i *number-of-jobs*) 
(print (elt (first *schedule*) i) out-stream) ) 
(print ‘Algorithm-Compute-Time out-stream) 
(print (- *algorithm-stop-time* *algorithm-start-time*) out-stream) 
(print ‘Meta-Task-Completion-Time out-stream) 
(print (get-completion-time) out-stream) 
(close out-stream} ) 


(defun dependency-status (child-job jobs-to-schedule) 
(let* ((status T) 
(parent-Jjob nil)) 
(dotimes (i (length *jobs*) status) 
(if (AND (= (elt (elt *dependency-matrix* i) child-job) 
1) 
(memberofp i jobs-to-schedule) ) 
(setf status NIL))))) 


(defun get-completion-time () 
(let* ((max-time 0Q)) 
(dotimes (i *number-of-jobs* max-time) 
(if (> (fourth (elt (first *schedule*) i)) max-time) 
(set£ max-time (fourth (elt (first *schedule*) i1))))))) 
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;iz27 Simulated annealing algorithm 


(defun set-job-list () 
(let* ((job-list nil)) 
(dotimes (i *number-of-jobs* job-list) 
(setf job-list 
(append job-list (list i)))))) 


(defun set-machine-list () 
(let* ((machine-list nil)) 
(dotimes (i *number-of-machines* machine-list) 
(setf machine-list 
(append machine-list (list i)))))) 


(defun get-time-to-compute (job machine) 
(elt (elt *etc-matrix* job) machine) ) 


(defun astarsearch (Start input-file output-file) 
(let* ((schedule nil)) 

(read-input-file input-file) 

(set£ *algorithm-start-time* (get-universal-time) ) 

(setf schedule (astarsearch2 

(list (list (evalfn start) start)) nil)) 

(set£ *algorithm-stop-time* (get-universal-time) ) 
(write-output-file output-file schedule) )) 


(defun astarsearch2 (agenda usedagenda) 
(when agenda 


(LE 

(goalreached (second (first agenda) )) 
(last (reverse (rest (first agenda) ))) 
({astarsearch2 

(merge 

‘list 

(remove-bad-agitems 

(mapcar 
#‘ (lambda 


(newstate) 
(let ((newpath 
(cons newstate (rest (elt agenda (get-state- 
to-visit agenda)))))) 
(cons (+ (evalfn newstate) (costfn newpath) ) 
newpath) ) ) 
(successors (second (elt agenda (get-state-to-visit 
agenda) )))) 
(append agenda usedagenda) ) 
(rest agenda) 
#’firstlessp) 
(cons (first agenda) usedagenda))))) 


(defun remove-bad-agitems (agitems oldagitems) 
(sort 
(remove-if 
#’ (lambda (agitem) 
(some #’ (lambda (oldagitem) 
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(and (>= (first agitem) (first oldagitem) ) 
(equalstate (second agitem) (second oldagitem) ))) 
oldagitems) ) 
agitems) F 
#’firstlessp) ) 


(defun firstlessp (agiteml agitem2) 
(< (first agiteml) (first agitem2))) 


(defun successors (state) 
(let* ((immediate-successors (Successors2 state) ) 
(complete-successors nil)) 
(dotimes (i (length immediate-successors) complete-successors) 
(set£ complete-successors 
(append complete-successors 
(list (append (elt immediate-successors 1) 

state))))))) 


(defun successors2 (state) 
(let* ((successors nil) 
(jobs-to-schedule (get-jobs-to-schedule state) ) 
(machine-stop-times (get-machine-stop-times state) )) 
(dotimes (m (length *machines*) successors) 
(dotimes (j (length jobs-to-schedule) nil) 
(if (dependency-status (elt jobs-to-schedule j) 
jobs-to-schedule) 
(setf successors 
(cons 
(list (list (elt jobs-to-schedule j) 
(elt *machines* m) 
(elt machine-stop-times m) 
(+ (elt machine-stop-times m) 
(get-time-to-compute (elt jobs-to-schedule 
j) m)))) 
successors))))))) 


(defun goalreached (statelist) 
(if (equalp (length (get-jobs-to-schedule statelist)) 0) t nil)) 


(defun equalstate (statel state2) 
(equal statel state2)) 


(defun costfn (statelist) 
(let* ((cost 0) 
(machine-cost (get-individual-machine-cost (first 
statelist)))) 
(dotimes (n (length machine-cost) cost) 
(if (> (elt machine-cost n) cost) 
(set£ cost (elt machine-cost n)))))) 


(defun evalfn (statelist) 
(if (= (length (get-jobs-to-schedule statelist)) 0) 0 
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‘@aea (>= (length (get-jobs-to-schedule statelist) ) (length 
*machines*) ) 
(/ (get-unscheduled-job-time statelist) (length *machines*) ) 
(/ (get-unscheduled-job-time statelist) (length (get-jobs-to- 
schedule statelist)))))) 


. (defun get-jobs-to-schedule (state) 
(let* ((jobs-to-schedule *jobs*)) 
(dotimes (n (length state) jobs-to-schedule) 
(if (memberofp (first (elt state n)) jobs-to-schedule) 
(set£ jobs-to-schedule 
(remove (first (elt state n)) jobs-to-schedule)))))) 


(defun get-average-job-time () 
(let* ((average-job-time (init-job-list))) 
(dotimes (i (length *jobs*) average-job-time) 
(dotimes (j (length *machines*) ) 
(setf (elt average-job-time i) 
(+ (elt average-job-time i) 
(get-time-to-compute i j)))) 
(setf (elt average-job-time 1) 
(/ (elt average-job-time i) (length *machines*)))))) 


{defun get-unscheduled-job-time (statelist) 
(let* ((time-left 0) 
(jobs-remaining (get-jobs-to-schedule statelist))) 
(dotimes (n (length jobs-remaining) time-left) 
(setf£ time-left 
(+ time-left (elt (get-average-job-time) (elt jobs-remaining 
n))))))) 


(defun get-individual-machine-cost (state) 
(let* ((machine-cost (init-machine-list) )) 
(dotimes (n (length state) machine-cost) 
(let* ((machine-to-update (second (elt state n))) 
(job-start-time (third (elt state n))) 
(job-finish-time (fourth (elt state n))) 
(job-scheduled (first (elt state n)))) 
(setf (elt machine-cost machine-to-update) 
(+ (elt machine-cost machine-to-update) 
(get-time-to-compute (first (elt state n)) 
(second (elt state n))))) 
(if (> job-finish-time (elt *deadline-matrix* job-scheduled) ) 
(set£ (elt machine-cost machine-to-update) 
(+ (elt machine-cost machine-to-update) | 
(- job-finish-time (elt *deadline-matrix* job- 
scheduled))))))))) 


(defun get-machine-stop-times (state) 
(let* ((machine-stop-times (init-machine-list))) 
(dotimes (i (length state) machine-stop-times) 
(let* ((machine (second (elt state i)))) 
(set£ (elt machine-stop-times machine) 
(fourth (elt state i))))))) 
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(defun init-job-list () 
(let* ((average-job-time nil)) 
(dotimes (i (length *jobs*) average-job-time) : 
(set£ average-job-time 
(append ‘(0) average-job-time))))) 


(defun init-machine-list () 
(let* ((machine-cost nil)) 
(dotimes (i (length *machines*) machine-cost) 
(setf£ machine-cost 
(append ‘(0) machine-cost))))) 


(defun notmemberofp (item list) 
(let* ((memberstatus t)) 
(dotimes (n (length list) memberstatus) 
(if (equalp item (elt list n)) 
(setf memberstatus nil))))) 


(defun memberofp (item list) 
(let* ((memberstatus nil)) 
(dotimes (n (length list) memberstatus) 
(if (equalp item (elt list n)) 
(set£f memberstatus t))))) 


(defun read-input-file (filename) 
(setq in-stream (open filename :direction :input) ) 
(let* ((number-of-jobs (read in-stream) ) 
(number-of-machines (read in-stream) ) ) 

(setf£ *number-of-jobs* number-of-jobs) 
(set£ *number-of-machines* number-of-machines) 
(setf *jobs* (set-job-list) ) 
(set£ *machines* (set-machine-list) ) 
(read-etc-matrix in-stream) 
(read-dependency-matrix in-stream) 
(read-deadline-matrix in-stream) 
(close in-stream) )) 


(defun read-etc-matrix (in-stream) 
(let* ((individual-job-times nil) 
(etc-matrix nil)) 
(dotimes (i *number-of-jobs*) 
(setf individual-job-times nil) 
(dotimes (j *number-of-machines* ) 
(setf individual-job-times 
(append individual-job-times (list (read in-stream))))) 
(setf etc-matrix 
(append etc-matrix (list individual-job-times) ))) 
(set£ *etc-matrix* etc-matrix) )) 


(defun read-dependency-matrix (in-stream) 
{(let* ((individual-job-dependencies nil) 
(dependency-matrix nil) ) 
(dotimes (i *number-of-jobs*) 
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(setf individual-job-dependencies nil) 
(dotimes (j *number-of-jobs*) 
(setf individual-job-dependencies 
(append individual-job-dependencies (list (read in- 
stream) )))) 
(set£ dependency-matrix 
(append dependency-matrix (list individual-job- 
a dependencies) ))) 
(setf *dependency-matrix* dependency-matrix) ) ) 


(defun read-deadline-matrix (in-stream) 
(let* ((deadline-matrix nil) ) 
(dotimes (1 *number-of-jobs*) 
(setf£ deadline-matrix 
(append deadline-matrix (list (read in-stream))))) 
(setf *deadline-matrix* deadline-matrix) )) 


(defun write-output-file (filename schedule) 
(setq out-stream (open filename :direction :output) ) 
(print ‘’LISP-Simulated-Annealing-Search-Summary out-stream) 
(dotimes (1 *number-of-jobs*) 
(print (elt (first schedule) i) out-stream) ) 
(print ‘Algorithm-Compute-Time out-stream) 
(print (- *algorithm-stop-time* *algorithm-start-time*) out-stream) 
(close out-stream) ) 


(defun dependency-status (child-job jobs-to-schedule) 
(let* ((status T) 
(parent-job nil)) 
(dotimes (1 (length *jobs*) status) 
(if (AND (= (elt (elt *dependency-matrix* i) child-job) 
1) 
(memberofp 1 jobs-to-schedule) ) 
(setf£ status NIL))))) 


(defun get-state-to-visit (agenda) 
(let* ((original-cost-list nil) 
(denominator 0) 
(weighted-cost-list nil) 
(normalized~-cost-list nil) 
(probability-window-list ‘(0)) 
(state-to-visit nil) 
(sum 0)) 
(setf original-cost-list (get-original-cost-list agenda) ) 
(setf denominator (get-denominator original-cost-list) ) 
(setf weighted-cost-list (get-weighted-cost-list original-cost- 
list) ) 
(setf annealed-cost-list (get-annealed-cost-list weighted-cost- 
* list) ) 
(set= normalized-cost-list (get-normalized-cost-list annealed- 
cost-list denominator) ) 


p (setf probability-window-list (get-probability-window-list 
normalized-cost-list) ) 
(sett state-to-visit (get-state-to-visit2 probability-window- 


LASt)) 
(setf sum (get-sum normalized-cost-list) ) 
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(print ‘weighted) 

(print weighted-cost-list) 

(print ‘annealed) 

(print annealed-cost-list) : 
(print ‘normalized) 

(print normalized-cost-list) 
(print ‘sum) 

(print sum) 

(print ‘probability) 

(print probability-window-list) 
(print ‘agenda-length) 

(print (length agenda) ) 

(print ‘state-to-visit) 

(print state-to-visit) )) 


(defun get-original-cost-list (agenda) 
(let* ((original-cost-list nil)) 
(dotimes (i (length agenda) original-cost-list) 
(sett original-cost-list 
(append original-cost-list 
(list (first (elt agenda 1)))))))) 


(defun get-denominator (original-cost-list) 
(let* ((denominator 0)) 
(dotimes (i (length original-cost-list) denominator) 
(setf£ denominator 
(+ denominator (/ (first original-cost-list) (elt original-cost- 


das. 332) 


(defun get-weighted-cost-list (original-cost-list) 
(let* ((weighted-cost-list nil)) 
(dotimes (i (length original-cost-list) weighted-cost-list) 
(setf weighted-cost-list 
(append weighted-cost-list 
(list (/ (first original-cost-list) (elt original- 


cost-list i)))))))) 


(defun get-annealed-cost-list (weighted-cost-list) 
(let* ((annealed-cost-list nil) 
(annealing-factor (get-annealing-factor) ) ) 
(if (< annealing-factor 1) 
{(setf annealing-factor 1)) 
(print ‘’annealing-factor) 
(print annealing-factor) 
(dotimes (i (length weighted-cost-list) annealed-cost-list) 
(if (= 1 Q) | 
(setf annealed-cost-list (list (first weighted-cost-list) )) 
(set£ annealed-cost-list 
(append annealed-cost-list 
(list (/ (elt weighted-cost-list i) annealing- 


factor)))))))) 


(defun get-normalized-cost-list (annealed-cost-list denominator) 
(let* ((normalized-cost-list nil)) 
(dotimes (i (length annealed-cost-list) normalized-cost-list) 
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(setf normalized-cost-list 
(append normalized-cost-list 
(list (/ (elt annealed-cost~list i} denominator))))))) 


(defun get-probability-window-list (normalized-cost-list) 
(let* ((probability-window-list /’(0))) 
(dotimes (i (length normalized-cost-list) probability-window-list) 
(if (= i 0) 
(setf probability-window-list 
(append probability-window-list 
(list (first normalized-cost-list) ))) 
(setf£ probability-window-list 
(append probability-window-list 
(list (+ (elt normalized-cost-list 1) 
(elt probability-window-list i))))))))) 


(defun get-state-to-visit2 (probability-window-list) 
(let* ((random-selector (random (float 1))) 
(state-to-visit nil)) 
(print random-selector) 
(dotimes (i (- (length probability-window-list) 1) state-to-visit) 
(eee 
(AND (if (> random-selector (elt probability-window-list i)) 
t) 
(if (< random-selector (elt probability-window-list (+ 
2 yy). El) 
(setf state-to-visit i))))) 


(defun get-sum (normalized-cost-list) 
(let* ((sum Q)) 
(dotimes (i (length normalized-cost-list) sum) 
(setf sum 
(+ sum (elt normalized-cost-list i)))))) 


(defun get-annealing-factor () 
(let* ({scaling-factor 1)) 
(* scaling-factor 
(- (get-universal-time) *algorithm-start-time*)))) 


77 


THIS PAGE INTENTIONALLY LEFT BLANK 


78 





LIST OF REFERENCES 


: [ALHUSAINI99] Ammar Alhusaini, Viktor Prasanna, and C.S. Raghavendra. A Unified 
Resource Scheduling Framework for heterogeneous Computing Environments. Proc. 8th 
EEE Heterogeneous Computing Workshop, Puerto Rico, April 1999. 


[ALI94] Hesham H. Ali, Theodore G. Lewis, and Hesham El-Rewini. Task Scheduling 
in Parallel and Distributed Systems. Prentice Hall, 1994. 


[ARMSTRONG98] Robert Armstrong, Debra Hensgen, and Taylor Kidd, The Relative 
Performance of Various Mapping Algorithms is Independent of Sizable Variances in 
Run-time Predictions, 7th Heterogenous Computing Workshop (HCW%8), Mar. 1998, 
pp. 79-87. 


[BHARADWAJ96] Veeravalli Bharadwaj, Debasish Ghose, Venkatarman Mani, and 
Thomas G. Robertazzi. Scheduling Divisible Loads in Parallel and Distributed Systems. 
IEEE Computer Society Press, 1996. 


[BHAT99] Prashanth B. Bhat, C.S. Raghavendra, and Viktor K. Prasanna. Efficient 
Collective Communication in Distributed Heterogeneous Systems. The 19th International 
Conference on Distributed Computing Systems (ICDCS), 1999. 


[(CORMEN97] Thomas H. Cormen, Charles E. Leiserson, and Ronald L. Rivest. 
Introduction to Algorithms. The MIT Press, 1997. 


[FREUND98] Richard F. Freund, Michael Gherrity, Stephen Ambrosius, Mark 
Campbell, Mike Halderman, Debra Hensgen, Elaine Keith, Taylor Kidd, Matt Kussow, 
John D. Lima, Francesca Mirabile, Lantz Moore, Brad Rust,and Howard Jay Siegel, 
Scheduling Resources in Multi-User, Heterogeneous, Computing Environments with 
SmartNet, 7th Heterogeneous Computing Workshop (HCW 98), Orlando, FL, Mar. 1998, 
pp. 184-199. 


[GALLIO0O] Doreen L. Galli. Distributed Operating Systems Concepts and Practice. 
Prentice Hall, 2000. 


[HURSON95] Ali Hurson, Krishna Kavi, and Behrooz Shirazi.” Scheduling and Load 
Balancing in Parallel and Distributed Systems. TEEE Computer Society Press, 1995. 


[KIDD96] Taylor Kidd, Debra Hensgen, Richard Freund, and Lantz Moore, SmartNet: A 
Scheduling Framework for Heterogeneous Computing, proceedings of the International 


3 symposium on Parallel Architectures, Algorithms, and Networks (I-SPAN%6), sponsored 
by the IEEE computer Society, Beijing, China, Jun. 1996, pp. 514-521. 


79 





[LUGER89] George F. Luger and William A. Stubblefield. Artificial Intelligence and 
the Design of Expert Systems. The Benjamin/Cummings Publishing Company, Inc., 1989. ; 


[MSHNO0] Management Systems for Heterogeneous Networks World Wide Web Page: 
http://www.mshn.org 











80 








rT 


BIBLIOGRAPHY 


: [ALI99] Shoukat Ali, A Comparative Study of Dynamic Mapping Heuristics for a Class 
of Independent Tasks onto Heterogeneous Computing Systems, Master’s Thesis, School 
of Electrical and Computer Engineering, Purdue University, 1999. 


[ARMSTRONG97] Robert Armstrong. Investigation of Effect of Different Run-Time 
Distributions on Smartnet Performance. Masters Thesis, Naval Postgraduate School, 
Monterey, CA, Sep. 1997. Debra Hensgen, advisor, Taylor Kidd, co-advisor. 


[BRAUN99] Tracy D. Braun, Howard Jay Siegel, Noah Beck, Ladislau L. Boloni, 
Muthucumaru Maheswaran, Albert J. Reuther, James P. Robertson, Mitchell D. Theys, 
Bin Yao, Debra Hensgen, and Richard F. Freund, A Comparison Study of Static Mapping 
Heuristics for a Class of Meta-tasks on Heterogeneous Computing Systems, Proceedings 
of the 8th IEEE Workshop on Heterogeneous Computing Systems (HCW 99), San Juan, 
Puerto Rico, Apr. 1999. 


[BRAUN98] Tracy D. Braun, Muthucumaru Maheswaran, Howard Jay Siegel, Noah 
Beck, Ladislau Boloni, Albert I. Reuther, James P. Robertson, Mitchell D. Theys, and Bin 
Yao, A Taxonomy for Describing Matching and Scheduling Heuristics for Mixed- 
Machine Heterogeneous Computing Systems, Workshop on Advances in Parallel and 
Distributed Systems (in proceedings of the IEEE Symposium on Reliable Distributed 
Systems, Oct. 1998, pp. 330-335), West Lafayette, IN, Oct. 1998. 


[KRUATRACHUE87] Boontee Kruatrachue, Static Task Scheduling and Grain Packing 
in Parallel Processing Systems, PhD Thesis, School of Electrical and Computer 
Engineering, University of Michigan, 1987. 


[KWOK99] Yu-Kwong Kwok, Anthony A. Maciejewski, Howard Jay Siegel, Arif 
Ghafoor, and Ishfaq Ahmad, Evaluation of A Semi-Static Approach to Mapping Dynamic 
Iterative Tasks onto Heterogeneous Computing Systems, 1999 International Symposium 
on Parallel Architectures, Algorithms, and Networks (I-SPAN %9), pp. 204-209, 
Fremantle, Australia, June 1999. 


[MAHESWARAN98] Muthucumaru Maheswaran, Tracy Braun, and Howard Jay Siegel, 
High-Performance Mixed-Machine Heterogeneous Computing, 6th Euromicro Workshop 
on Parallel and Distributed Processing, Madrid, Spain, Jan. 1998, pp. 3-9. 


[MAHESWARAN98] Muthucumaru Maheswaran and Howard Jay Siegel, A Dynamic 
Matching and Scheduling Algorithm for Heterogeneous Computing Systems, 7th 
‘ Heterogeneous Computing Workshop (HCW 98), Orlando, FL, Mar. 1998, pp. 57-69. 


[MAHESWARAN 99] Muthucumaru Maheswaran, Shoukat Ali, Howard Jay Siegel, 


Debra Hensgen, and Richard F. Freund, Dynamic Mapping of a Class of Independent 
Tasks onto Heterogeneous Computing Systems, Journal of Parallel and Distributed 


81 





Computing, Special Issue on Software Support for Distributed Computing, 1999, 
accepted and scheduled to appear. | 


[MAHESWARAN99] Muthucumaru Maheswaran, Tracy D. Braun, and Howard Jay 
Siegel, Heterogeneous Distributed Computing, Encyclopedia of Electrical and 
Electronics Engineering, J. G. Webster, editor, John Wiley & Sons, New York, NY, 
1999, Vol 8, pp. 679-690. 


[MAHESWARAN99] Muthucumaru Maheswaran, Shoukat Ali, Howard Jay Siegel, 
Debra Hensgen, and Richard F. Freund, Dynamic Matching and Scheduling of a Class of 
Independent Tasks onto Heterogeneous Computing Systems, proceedings of the 8th 
Workshop on Heterogeneous Computing Systems (HCW 99), San Juan, Puerto Rico, 


Apr. 1999. 


[MITCHELL96] Melanie Mitchell. An Introduction to Genetic Algorithms. The MIT 
Press, 1996. 


[ROSEN92] Benjamin Rosen. Advances in Optimization and Parallel Computing. 
North-Holland, 1992. 


[WANG97] Lee Wang, Howard Jay Siegel, Vwani P. Roychowdhury, and Anthony A. 
Maciejewski, Task Matching and Scheduling in Heterogeneous Computing Environments 
Using a Genetic-Algorithm-Based Approach, Journal of Parallel and Distributed 
Computing, Special Issue on Parallel Evolutionary Computing, Vol. 47, No. 1, Nov. 
1997, pp. 8-22. 


82 





INITIAL DISTRIBUTION LIST 


Defense Technical Information Center .00.........cccccccccccccccccccccecceccceccececeececccccccccecccece, 
8725 John J. Kingman Road, Suite 0944 
Ft. Belvoir, VA 22060-6218 


PSY FOX TSA aos shea yaeecnass casuavensesichagates oouscotedebav nae taensasthcancewesieesbameeniete 
Naval Postgraduate School 

411 Dyer Road 

Monterey, CA 93943-5101 


PLOLeSSOr Neth ROW sicstecees coco ee ata tvosheSasnce lei casi soasSeaebaetiablataadlacasehsstavarstense, 
Computer Science Department, Code 32 

Naval Postgraduate School 

Monterey, CA 93943-5101] 


Wiehe 1D NIC MOI csc sorsa cs etcoi.ic5 5 iileSayatwau ess dieucuaitanse setece Gertie cestode occas ae Sicietnest 
3908 Knoll Ridge Drive 
Cedar Falls, Iowa 50613 


Chairmain, Computer Science Department ...........0....cccscscescccesscceesecsccceessecessnsceees 
Code CS 

Naval Postgraduate School 

Monterey, CA 93943-5101 


83 





